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ABSTRACT 


This  thesis  describes  a  procedure  for  computing  a  lower  confidence 
limit  on  the  reliability  of  a  quasi-coherent  complex  system  using  test  data 
on  its  components.  The  failure  times  of  the  components  are  assumed  to 
have  either  exponential  or  Weibull  distributions  with  unknown  parame¬ 
ters.  The  accuracy  of  this  procedure  is  evaluated  using  computer  simu¬ 
lation  for  various  system  structures  and  sets  of  parameter  values  for  the 
assumed  distributions.  This  thesis  is  an  extension  of  a  thesis  by  Kah  Chee 
Yee  in  that  it  uses  a  different  equation  for  the  estimate  of  the  shape  pa¬ 
rameter  in  the  Weibull  distribution  than  Yee  used,  and  it  evaluates  the 
procedure  for  a  larger  collection  of  system  structures. 
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J.  INTRODUCTION 


This  thesis  modifies  an  existing  approximate  lower  confidence  limit 
procedure  for  the  reliability  of  complex  systems  developed  by  Yee 
CRef.  73.  The  purpose  of  the  modification  is  to  make  the  procedure  more 
accurate  fijr  more  configurations  of  quasi-coherent  complex  systems  and 
easier  to  use  computationally.  A  system  is  defined  to  be  quasi-coherent 
if  an  increase  in  reliability  of  any  one  of  its  components  does  not  cause 
a  decrease  in  the  system  reliability.  The  components  of  a  quasi-coherent 
system  do  not  need  to  be  statistically  independent.  However,  throughout 
this  thesis  it  is  assumed  that  all  components  are  statistically  independent 
and  the  probability  distributions  of  their  failure  times  are  either  expo¬ 
nential  or  Weibull . 

The  approximate  system  reliability  lower  confidence  limit  procedure 
developed  by  Yee  is  quite  accurate  for  series  systems.  In  this  thesis  the 
procedure  for  estimating  the  shape  parameter,  P,  differs  from  the  maxi¬ 
mum  likelihood  method  used  by  Yee.  The  procedure  developed  by 
Varadan  CRef.  6]  is  used  to  estimate  the  shape  parameter,  in  the 
Weibull  distribution  in  this  thesis. 

In  addition  to  modifying  the  lower  confidence  interval  estimation 
procedure  developed  by  Yee,  more  complex  structures  are  examined  here 
than  in  Yec's  thesis.  The  computer  programs  developed  by  Yee  v/ere 
modified  to  examine  these  new  structures.  Also,  a  computer  program  was 
developed  that  can  be  used  to  compute  the  lower  confidence  limit  for  the 
reliability  of  a  complex  system  using  these  procedures.  The  system  is  de- 
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scribed  by  the  user  in  response  to  queries  by  the  program.  Also,  the  test 
data  set  :s  entered  by  the  user  in  response  to  queries. 
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II.  ESTIMATES  FOR  THE  SHAPE  PARAMETER  OF  A  WEIBULL 

DISTRIBUTION 

If  the  time  to  failure,  X,  of  a  device  has  a  Weibull  distribution  with 
scale  and  shape  parameters  Q  and  ^  respectively,  then  its  probability 
density  function  is  given  by 

g{x;8,^)=^{\iefpx^~^exp{-{xldf}  ,  x>0  (2.1) 

This  property  is  stated  more  briefly  by  the  phrase  X  is  WEI(  0,  P  ).  The 
cumulative  distribution  function,  CDF  for  X  is 

Gix;  ej)^\  -  exp{  -  {x/df  ]  ,  x  >  0  (2.2) 

and  the  survival  function  is 

G(x;0,^)=l  -  G{x;0J)  (2.3) 

It  is  well  known,  see  [Ref.  63,  that  the  random  variable  T=  In  Zhas  the 
extreme  value  distribution  with  density  function  f  and  CDF  F  given  by 

Jljf;  n,  a)  =  exp[  -  oo  <  y  <  oo 

~  oo  <  n  <  oo  (2"4) 
a  >  0 

F(y\  //,  <t)  =  1  -  exp{  (2.5) 

where  o  =  1//? ,  ju  =  In  6  and  In  is  the  natural  logarithm.  In  this  case  we 
write  Y  is  EV(  a,  fi ). 
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Suppose  Z,,Z2,...,Z„  are  independent  random  variables  with  a  common 
WEI(  6,^)  distribution  and  let  <  Xf^)  <  ...  <  X^^)  be  the  first  r  ordered 
set  of  the  original  set  of  n  variables-  Then  the  following  equations  are 

A  A 

solved  to  obtain  the  maximum  likelihood  estimates  (MLE)  d  and  (j  for  6 
and  P,  see  [Ref.  l]  : 

£  4  In  +  (n  -  In  Aj,, 

J,  4 +(«-'■<) 


JL 

P 


In  X( 


V) 


(2.6) 


1 


_ r _ 

4-) + (« -^'^4) 


(2.7) 


Closed  form  expressions  for  p  do  not  exist,  and  iterative  procedures  are 

A  A 

used  to  solve  for  P  and  6.  Computer  programs  are  readily  available  to 

A  A 

compute  P  and  6.  They  require  an  original  estimate  Po  to  start  an  iterative 
process,  and  if  it  is  not  chosen  carefully  the  iteration  will  not  yield  an 

A 

accurate  estimate  for  p.  The  MLE  estimate,  P,  is  biased.  Bain  and  Lee 
[Ref.  5]  have  an  excellent  treatment  of  the  Weibull  distribution  and 

A  A 

properties  of  the  MLE  estimates  6  and  p.  Their  discussion  includes  a  table 
of  factors  (page  200)  which  can  be  used  to  construct  nearly  unbiased  es¬ 
timates  for  P  . 

Balakrishnan  and  Varadan  [Ref.  6]  have  derived  a  method  for  esti¬ 
mating  p  and  6  which  does  not  require  computer  iterations.  Their  method 
will  be  used  in  this  thesis  to  modify  Yee's  existing  procedure  that  derives 
lower  confidence  limits  for  the  reliability  of  complex  systems  that  have 
some  Components  whose  failure  times  have  a  Weibull  distribution.  A 
summary/  of  their  results  follows. 
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Let  5^(r+i)iL'(,^.2),..M  .  .  '--e  a  sequence  of  the  order  statistics  for  a  ran¬ 

dom  sample  of  size  n  variables  each  with  density  function  given  by  (2.4). 


Then 


^(r+l)  ^  Y ;r+2)  ^  ^  (2.8; 

If  we  let  s  —  fi)f a,  then  the  likelihood  function  based  on  the 

censored  sample  is : 

i Cf(4+i))T  p.9: 

where, 

A  =  n~  r  —  s 
Fix)  s  1  -  exp(  -(Y) 

Fix)  s  1-F(x) 
fix)  s  exp(  -e^) 

The  likelihood  equations  for  n  and  a  are  : 


^InL 

Sn 


=-^c 


1  j-  fi^{n-s))  ^ 


n—5  ft 


]  =  0  (2.10) 


5\nL 

do 


(2.U) 


•^"^(w-5)  pfY  )  ^iO  f(X  )  ^  ~ 
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Equations  (2.10)  and  (2.11)  do  not  admit  explicit  solutions.  Let,  denote 
the  approximated  likelihood  function.  Using  a  Taylor  series  expansion  for 
In  L  to  obtain  In  L’‘,  the  two  equations  can  be  approximated  by  : 


5\nL  5  Ini'" 

- - ^ - 

Oli  on 


(2.12) 


n—s 


5(1  -  (a,  -  ^,Z(,-))]  =  0 


^InL  ^InL^ 
<5(7  da 


(2.13) 


•^^(rt-s)(^  “  ^n-s  +  Pn-s'^n-s))  ~  “  0* 

See  Balakrishnan  and  Varadan  [Ref.  6  p-147] 

Solving  (2.12)  and  (2.13)  for  ju  and  a  gives  their  approximate  maxi¬ 
mum  likelihood  estimates  as  follows  : 


n  =  B-aC 

(2.14) 

a  =  {-D  +  iD^  +  4AE)^^^]I2A 

(2.15) 

where  , 

_  i 

“  {n+\) 

(2.15) 

h'i  =  1-P/ 

(2.17) 

a,-  =  l-h  In  qi{\  -  ln(  -  In  ql)] 

(2.18) 
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Pi  = 


(2.19) 


S 


^r+l 

Pr+l 


In  ^^+1} 


(2.20) 


^y+l 

Pr+\ 


In  -  ln(  -  In  ^^^.i)} + 


(2.21) 


_fr+l 

2  “ 
Pr+l 


(ln^;-+l) 


2 


ln(  -  In  ^^+1) 


A^n—r  —  s  (2.22) 

B  =  {fd  ]/m  (2.23) 

C={r)'-5(1 ’-a„__j)+  S  «/  ]/w  (2.24) 

/=r+l 

m  =  +sPfj-s  "I"  ^  Pi  (2.25) 

/=/'+l 

O  =  '■’X^ir+D  -.B)  -  XI  -  -B)  -S)  (2.26) 

■E=  '•.5(7(,+i)  -Bf  +  JjS„_,(y(„_s)  -  Xj+i  W^(0  -  P-2’) 


The  approximate  MLE  estimates  of  fi  and  a  as  defined  in  equations  (2.14) 
and  (2.15)  are  biased  estimators.  Balakrishnan  and  Varadan 
LRef.6  p  -  149]  provide  a  table  of  constants  that  are  called  SBIAS(n,r,s). 
They  are  a  function  of  the  number  n  of  test  items  placed  on  test,  the 
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number  r  where  the  observations  began  and  the  number  s  of  successive 
components  observed.  In  this  thesis  the  parameter  r  will  always  be  zero, 
meaning  that  observations  always  start  from  the  smallest  failure  time. 
Using  their  table  an  approximate  unbiased  estimate  d*'  for  cr  is  (See  Ap¬ 
pendix  A) 

^  “  ~SBIAS{n,0,s) 

^  A 

The  inverse  of  a  is  jS  , which  will  be  a  biased  estimator  for  /?  .  Bain 
[Ref.  5  p-220]  provides  a  table  of  constants  B{n),  which  depend  on  the 

A  A 

number  of  test  items,  n,  such  that  =  j?  {5(«)}  is  nearly  unbiased  for  p. 
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III.  DESCRIPTION  OF  THE  LOWER  CONFIDENCE  LIMIT 

PROCEDURE 


The  lower  confidence  limit  procedure  developed  in  this  thesis  is  a 
modification  of  the  procedure  in  the  thesis  written  by  Yee  ,  [Ref.  7]  . 
The  procedure  in  this  thesis  uses  different  estimators  for  the  parameters 
P  and  6  in  the  Weibull  distribution.  The  method  of  estimation  presented 
by  Mann  and  others  [Ref.  l]  and  used  by  Yee  employs  the  maximum 
likelihood  estimates  for  p  and  6  which  require  computer  iteration  and  a 
reasonably  good  guess  for  an  initial  value  of  P  to  perform  the  iteration. 
The  estimators  developed  by  Balakrishnan  and  Varadan  [Ref.  6]  is  used 
in  this  thesis.  They  provide  an  alternative  estimation  procedure  which 
does  not  require  computer  iteration  methods  to  compute  the  estimates. 

In  this  thesis  we  consider  systems  that  are  made  up  of  k  independent 
component  subsystems.  Systems  undergo  missions  of  some  duration,  say 
t.  During  a  mission  the  system  components  are  subject  to  periods  of  ac¬ 
tivity  and  inactivity.  During  active  periods  component  i  is  subject  to 
failure,  with  a  hazard  rate  4(/).  During  inactive  periods  the  hazard  rates 
are  all  zero.  A  component  successfully  completes  a  mission  if  its  total 
operating  time  in  [O,  r]  exceeds  some  design  parameter  t/r).  A  system 
completes  its  mission  if  sufficiently  many  of  its  subsystems  do;  the  system 
reliability  is,  as  usual,  a  function  of  the  structure  of  the  system. 

A.  SERIES  SYSTEMS 

1.  Inten  al  Estimation  Procedure  for  Exponential  Failure  Times 
Suppose  a  series  system  has  k  components  whose  failure  times  are 
statistically  independent.  Suppose  the  distribution  of  the  failure  time  of 
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component  i  is  exponential  with  failure  rate  Then  the  system  reliability 
can  be  written  as  a  function  of  Xi  and  i  =  l,2,...,k  as  follows  : 


k 

i?_j(r)  =  exp(  -S  Xiti) 

i=l 


(3.1) 


where  ti  =  ///)  is  the  time  component  i  operates  when  the  system  operates 
for  time  t.  Using  the  relationship  ri  =  X-JX^,  for  i=  where  2^  is 

the  failure  rate  of  any  one  of  the  k  components,  equation  (3.1)  becomes 

Rlt)  =  exp(  -  X^J.  Vjti )  (3.2) 

/=1 

A 

If  the  values  of  the  r,  are  known  and  if  is  an  upper  100(1  —  a)  % 
confidence  limit  for  X„,  the  corresponding  approximate  lower  confidence 
limit  for  R,{t)  would  be  : 

A  A  /: 

^/^LCa)  =  exp(  -  />.^,(;(a).5j  '•/?/ )  (3*3) 


The  equation  for  2„  depends  on  the  plan  for  testing  the  components. 
The  following  case  is  considered  in  this  thesis  : 

If  rii  items  of  component  type  i  are  tested  until  fi  failures  occur,  i 
=  l,2,...,k  and  if  ^(d,  denote  these  ordered /  failure  times  then 


2 

^  ^a,2F 

^'m,U{a)  ~  I 

2  S  nT, 

l=l 


(3.4) 


where  T,  denotes  the  total  test  time  accumulated  on  all  «,•  items  of  t5T3e 
i;  i.e,  7)  =  (n, ->9%  +  F  =  H  and  xhp  is  the  100(1 -a)^/2 

J=\  /=! 
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percentile  point  of  a  Chi-square  distribution  with  2F  degrees  of  freedom. 
See  Bain  and  Engelhardt  [Ref.  3  ]. 

Values  of  the  r,  are  assumed  to  be  unknown  in  this  thesis.  They 
will  be  estimated  by  r, ,  a  nearly  unbiased  estimator  for  r, ,  defined  by 


^'nt 

A 

where  ?.i  =  (fi—\)ITi  which  is  an  unbiased  estimator  for  when  fi>  1 

_  _  A  A  A 

[Ref.  7  p.  39 J.  If  were  unbiased  for  l/2„  and  if  and  2,  are  inde- 

pendent,  then  r,-  would  be  an  unbiased  estimator  for  r^.  Replacing  with 
equation  (3.5)  will  yield  an  unbiased  estimator  r,  for  ?/. 
Multiplying  by  this  constant  —1)  is  nullified  by  a  cancellation  with 
the  same  constant  in  the  final  equation  for  the  system  reliability  lower 
confidence  limit,  so  equation  (3.5)  is  used  to  estimate  r,.  Using  the  esti¬ 
mator  %  for  r„  equation  (3.4)  becomes 


2 

A  Xol,2F 

^'m,U{a)  ~  k  ^ 

21  riT: 
i=l  ‘  ‘ 


(3.6) 


It  is  important  to  note  that  the  index  m  denotes  the  component  for  which 

A 

/,  =  —  is  the  largest  among  all  the  components  in  the  system.  The 

corresponding  equation  for  the  100(1  —  a)  %  lower  confidence  limit  on 
the  reliability  of  the  series  system  becomes 

4(^)L(a)  =  exp(  -  4,{/(a)  )  (3-7) 
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2.  Interval  Estimation  Procedure  for  Weibull  Failure  Times 
Suppose  a  series  system  has  k  statistically  independent  components 
and  suppose  the  failure  times  of  all  k  components  have  Weibull  distrib¬ 
utions.  Let  the  time  to  failure,  ^  ,  of  components  of  type  i  have  density 
function 


m  expt  -  1  ,  !;  >  0  (3.8) 

Then  the  system  reliability  i?//)  can  be  vi^ritten  as 

=  (3.9) 

J=1 


where  Af  =  xf' ,  is  any  one  of  the  Xf  ,  i=  l,2,...,k  and  ri  —  .  If  the 

values  are  known,  then  Xf>  will  have  an  exponential  distribution  with 
constant  failure  rate  and  the  procedures  described  in  Section  1  can  be 

A 

used  to  obtain  with  T,  defined  by 

f 

=  («,-y34) +;i,  4  P-'O) 

Suppose  Pi  is  unknown  and  are  the  ordered  failure  times 

under  failure  truncated  testing  for  component  i  in  the  system.  The  maxi- 

A  A 

mum  likelihood  estimates  pi  for  Pi  and  2,  for  2,  are  the  solutions  for  Pi 
and  Xi  in  equations  ; 
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X 

Pi 


ill  “* 


(3.11) 


Af  =  - - -  (3.12) 

JjA;§,  +  («,.-y34) 

where  under  failure  truncation  testing.  See  Mann  and  others 

[Ref.  1  p.  189-191]  . 

Equations  (3.11)  and  (3.12),  can  not  be  solved  in  closed  form. 
They  may  be  solved  using  iterative  computer  methods,  but  in  this  thesis 
another  method  is  used  which  does  not  need  computer  iterations  to  ob- 

A  A 

tain  15 1  and  2,  .  This  method  was  described  in  Chapter  II.  It  is  an  ap¬ 
proximation  of  the  maximum  likelihood  estimate.  Let, 

Tij  =  jq  ,  z=  l,2,...,k  y=  1,2,...,  (3.13) 


In  this  thesis  the  distribution  of  7]y  is  approximated  by  the  exponential 

A 

distribution  with  failure  rate  =  Xf. 

The  procedures  for  obtaining  the  lower  confidence  limit  for  system 
reliability  are  similar  to  those  in  Section  A.  Define 


where, 
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fi  _  fi 
T:  «/ 

j=\  J 


,  i~  1,2,. ..,/c 


(3.15) 


2^  =  max  . 
i 


Then  equation  (3.14)  becomes 


A  'iX/  '^m  ^ 

Jm 


(3.16) 


The  approximate  100(1  —  a)  %  upper  confidence  limit  for  is 


^  A 

2  S  r,r, 

i=i 


(3.17) 


where. 


=  1  f:  . 
i=V‘ 


(3.18) 


The  corresponding  approximate  100(1  —  a)  %  lower  confidence  limit  for 
the  reliability  7^(0  of  the  series  system  is  given  by  : 


(3.19) 


The  accuracies  of  these  approximate  confidence  interval  procedures  were 
evaluated  by  using  computer  simulations  which  are  described  in  the  next 
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chapter.  During  this  evaluation  process,  the  degrees  of  freedom  in  the 
expression  in  equation  (3.17),  was  increased  and  decreased  from  the 
defined  values  of  P  given  by  this  equation.  The  purpose  of  these  mod¬ 
ifications  v;as  to  find  more  accurate  lower  confidence  limit  procedures. 
The  specific  increases  and  decreases  are  described  in  chapter  IV.  The  re¬ 
sults  show  that  for  some  cases  the  procedure  with  modified  degrees  of 
freedom  is  more  accurate. 

B.  PARALLEL  SYSTEMS 

An  active  parallel  one-out-of-k  system  is  defined  as  a  system  consisting 
of  k  subsystems  such  that  system  failure  occurs  when  and  only  when  all 
k  subsystems  fail.  Equivalently,  the  system  is  successful  when  at  least  one 
of  its  subsystems  is  successful.  Such  a  parallel  system  is  said  to  be  an  ac¬ 
tive  redundant  system  of  order  k.  It  can  be  shown  that  all  of  the  prop¬ 
erties  given  in  the  preceding  section  for  serial  systems  can  be  dualized  to 
the  corresponding  properties  for  parallel  systems  by  replacing  any  event 
by  its  complementary  event. 

Let  the  parallel  system  be  made  up  of  k  independently  operating 
components,  each  with  an  exponentially  distributed  failure  time  and  fail¬ 
ure  rate ,i=  \2,...,k.  Then  the  system  reliability,  Rif),  can  be  written 
as  a  function  of  ,  i-  l,2,...,k  as  follows  : 

-  I  -  n  r  1  -  exp(  -  Xiti )  ]  (3.20) 

Using  equations  that  have  been  derived  in  the  preceding  section;  i.e., 
in  the  equations  (3.5),(3.6),(3.14),(3.15),(3.16)  and  (3.17),  the  correspond¬ 
ing  equations  for  the  approximate  100(1  —  a)  %  lower  confidence  limit 
for  the  reliability  of  a  parallel  system  is 
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for  exponential  failure  times. 

For  Weibull  failure  times  the  approximate  100(1  ~  a)  %  lower  confi¬ 
dence  limit  for  the  reliability  of  a  parallel  system  is 


(3.22) 


C.  SERIES-PARALLEL  SYSTEMS 

Series-parallel  or  parallel-series  systems  come  in  many  varieties.  One 
example  of  an  active  series-parallel  system  is  defined  as  a  system  that  is 
comprised  of  k  subsystems  in  series  connected  to  d  subsystems  in  parallel 
or  vice-versa.  The  reliability  block  diagram  of  an  example  system  can  be 
seen  in  Figure  3.1. 


Fig.3.1  Block  Diagram  Of  Series-parallel  System 


The  system  reliability  R/j)  of  this  series-parallel  system  can  be  written 
Rlt)  =  [  exp(  -  1  X;ti )] C 1  -  n  [1  -  exp(  -  ;.///)]]  (3.23) 

,•=1  t=i 

if  the  failure  time  of  each  component  is  exponential,  and 

R/t)  =  [  exp(  -  I  Xl'if' )][  1  -  n  (1  -  exp(  -  P.f'tf' )]]  (3,24) 

/=!  /=! 

if  the  failure  time  of  each  component  has  a  Weibuli  distribution. 

Using  the  same  approach  as  those  used  for  series  systems  and  parallel 
systems  for  estimating  the  parameters,  the  corresponding  equations  for 
an  approximate  100(1  —  a)  %  lower  confidence  limits  for  the  reliability 
of  the  series-parallel  system  are 

t  ■”  0^/  )^*  (3.25) 

■U  -T^n  ^  [1  ~  exp(  -  Piti)]] 


and 


4(^)L(a)  =  n  exp( 


^'m,  U(a)  J  -I  • 


(3.26) 


IV.  COMPUTER  SIMULATIONS 


TEST  PLAN  :  TESTING  COMPONENTS  UNTIL  fi  FAIL  (RETP) 

RETP  is  a  program  written  in  FORTRAN,  on  the  Amdahl  mainframe 
computer.  It  performs  the  computer  simulations  that  assess  the  accuracy 
of  the  lower  confidence  limit  procedure  for  system  reliability.  A  doc¬ 
umentation  of  this  program  and  its  associated  subroutines  is  included  in 
Appendix  C. 

The  program  accepts  input  parameters  via  an  input  tile  INI.DAT. 
For  each  replication,  it  generates  the  failure  times  for  all  the  component 
items  included  in  the  test  plan  using  a  uniform  random  number  generat¬ 
ing  subroutine  LRNDPC.  The  program  determines  the  total  test  time 
accumulated  for  2ach  component  in  the  system  and  computes  the  esti¬ 
mates  of  the  key  parameters  and  the  consequent  lower  confidence  limit 
for  system  reliability  for  that  replication.  This  is  done  for  various  system 
configurations;  i.e.,  series  systems  (Fig.  4.1),  parallel  systems  (Fig.  4.2), 
parallel-series  systems  (Fig.  4.3)  and  a  more  general  series-parallel  system 
(Fig.  4.4).  For  each  specific  system  configuration  and  set  of  input  pa¬ 
rameters,  the  process  is  repeated  1000  times.  When  all  replications  are 
done,  the  routine  EVAL  processes  the  lower  confidence  limit  estimates 
from  all  1000  replications  and  determines  the  two  measures  of  accuracy 
for  the  run,  namely  RSLOW  and  LEVEL. 

RSLOW  is  the  100(1  —  a)  percentile  of  the  ordered  set  of  lower  con¬ 
fidence  limits  from  the  1000  replications  computed  in  a  run.  The  true 
reliability  of  the  system  is  RS.  The  closer  RSLOW  is  to  RS,  the  greater 
the  accuracy  of  the  procedure  under  evaluation  in  the  run.  If  the  proce- 


18 


dure  is  exact,  RSLOW  will  be  equivalent  to  RS.  To  be  conservative, 
RSLOW  should  be  lower  than  RS. 

LEVEL  measures  the  proportion  of  1000  lower  confidence  limits, 
from  a  run  with  1000  replications,  which  are  lower  than  the  true  system 
reliability  RS.  The  closer  LEVEL  is  to  the  specified  confidence  ?  wel  for 
the  procedure,  1  -  a,  the  better  the  procedure.  Values  of  LEVEL  greater 
than  (I  —  a)  reflect  an  under-estimation  of  RS  which  is  conservative. 
Values  of  LEVEL  less  than  1  —  a  signal  an  over-estimation  of  RS  which 
may  be  undesirable. 


Fig.  4.1  Block  Diagram  Of  Series  System 


Fig.  4.2  Block  Diagram  Of  Parallel  System 


Fig.  4.3  Block  Diagram  Of  Parallel-series  System 
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Simulation  runs  are  performed  using  RETP  for  all  combinations  of 
failure  time  distributions  and  levels  of  key  input  parameters  listed  below. 
1.  System 

a.  8  Exponential  (Exp)  components  in  Series  (Case  1) 

b.  8  Weibull  (Wei)  components  in  Series  (Case  2) 

c.  4  Exp  and  4  Wei  (Mixed)  components  in  Series  (Case  3) 

Results  from  these  simulation  runs  are  put  into  tabular  form. 
By  observing  the  most  accurate  result  for  different  degrees  of 
freedom,  an  approximate  equation  for  obtaining  the  Degrees  of 
Freedom  can  be  derived,  see  Appendix  B.  Using  this  equation  we 
then  performed  simulation  runs  for  other  systems  for  RETP. 

d.  8  Exponential  components  in  Parallel  (Case  4) 

e.  8  Weibull  components  in  Parallel  (Case  5) 

f.  4  Exp  and  4  Wei  (Mixed)  components  in  Parallel  (Case  6) 

g.  8  Exponential  components  in  Parallel-series  (Case  7) 

h.  8  Weibull  components  in  Parallel-series  (Case  8) 
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i.  4  Exp  and  4  Wei  (Mixed)  components  in  Series-parallel  (Case  9) 

j.  5  Exp  and  10  Wei  (Mixed)  components  in  Series-parallel  (Case  10) 

2.  True  System  Reliability  (RS) 

a.  Hi  (greater  than  0.9)  (Type  A) 

b.  Lo  (greater  than  0.8)  (Type  B) 

See  pages  9-11  for  description  of  the  values  of  X;  and  for  each  of 
the  above  10  cases. 

3.  Level  of  Significance  (  a  ) 

a.  0.1 

b.  0.2 

4.  Degrees  of  Freedom  (DF)  for  the  statistic  as  a  function  of  the  total 
number  of  failed  test  components  (NFC)  and  total  number  of  system 
components  (NCOMP),  (For  cases  1,2  and  3  only). 

a.  DF  =  2NFC 

b.  DF  =  2(NFC  +  NCOMP) 

c.  DF  =  2(NFC  -  NCOMP) 

d.  DF  =  2NFC  -  NCOMP 

For  case  4  up  to  10  the  equation  derived  in  Appendix  B  was 
used.  This  equation  is 

e.  DF  =  2NFC  +  0.5(NC/NF) 

5.  Test  Plan  for  each  component. 

a.  Test  10  until  10  failures 

b.  Test  15  until  15  failures 

c.  Test  15  until  1 1  failures 

d.  Test  15  until  7  failures 

e.  Test  15  until  3  failures 

For  the  8  exponential  components  in  series  in  Case  1,  the  mission  time 
for  each  of  the  components  was  chosen  to  be  10  hours.  The  chosen  values 
of  the  scale  parameters,  A, ,  are  different  depending  on  whether  the  system 
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is  highly  reliable  (Type  A  system)  or  one  with  a  lower  reliability  (Type  B 
system).  The  ratios  between  the  largest  and  the  smallest  failure  rate  were 
chosen  to  be  8  and  4.5  respectively  for  Type  A  and  Type  B  systems. 

For  the  8  Weibull  components  in  series  in  Case  2,  the  mission  time 
for  each  of  the  components  was  chosen  to  be  10  hours.  In  order  to  obtain 
a  highly  reliable  system  (Type  A)  or  lower  reliability  system  (Type  B)  the 
scale  parameters,  2/,  were  chosen  differently.  The  ratios  between  the  larg¬ 
est  and  the  smallest  failure  rate  were  chosen  to  be  8  for  both  system  types. 
The  shape  parameter,  was  chosen  to  be  2  for  all  subcases.  The  program 
will  accomodate  any  value  greater  than  zero  for  the  shape  parameter. 

In  Case  3,  which  is  a  mixture  of  exponential  and  Weibull  compo¬ 
nents,  the  mission  time  for  each  of  the  components  was  chosen  to  be 
equal  to  10  hours.  The  scale  parameter  for  each  component  type  was 
chosen  so  that  the  ratios  between  the  largest  and  the  smallest  failure  rate 
were  4  for  both  component  types.  The  shape  parameter  for  each  of  the 
Weibull  components  is  chosen  to  be  2  for  all  cases. 

In  Case  4(8  exponential  components  in  parallel ),  the  mission  time 
for  each  of  the  components  was  chosen  to  be  10  hours.  The  chosen  values 
of  the  scale  parameters,  2^ ,  were  different  depending  on  whether  the  sys¬ 
tem  is  highly  reliable  (Type  A  system)  or  one  with  a  lower  reliability 
(Type  B  system).  The  ratios  between  the  largest  and  the  smalllest  failure 
rate  were  chosen  to  be  1.525  and  2.575  respectively  for  Type  A  and  Type 
B  systems. 

For  the  8  Weibull  components  in  parallel  in  Case  5,  no  change  was 
made  for  the  mission  time  for  each  of  the  components  as  described  in  the 
previous  paragraph.  The  ratios  between  the  largest  and  the  smallest  failure 
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rate  were  chosen  to  be  1.28  for  Type  A  systems  and  1.63  for  Type  B  sys¬ 
tems.  The  shape  parameter,  /?„  was  chosen  to  be  2  for  all  subcases. 

For  a  mixture  of  exponential  and  Weibull  components  in  parallel, 
the  mission  time  is  still  the  same  as  before  which  is  10  hours.  The  values 
of  the  scale  parameters,  for  the  exponential  components  were  chosen 
to  be  equal  to  those  for  the  Weibull  components.  The  ratios  between  the 
largest  and  the  smallest  failure  rates  were  chosen  to  be  l.j  and  2.2  for 
Type  A  and  Type  B  systems  respectively.  The  shape  parameter,  j?/,  for  a 
Weibull  component  was  chosen  to  be  2  for  all  subcases. 

In  Case  7,  the  8  exponential  components  in  parallel-series  (Figure 
4.3),  the  mission  time  for  each  of  the  components  was  chosen  to  be  10 
hours.  The  chosen  values  of  the  scale  parameters,  A,  were  different  de¬ 
pending  on  whether  the  system  is  highly  reliable  (Type  A  system)  or  one 
with  a  lower  reliability  (Type  B  system).  The  ratios  between  the  largest 
and  the  smallest  failure  rate  were  chosen  to  be  8  for  both  Type  A  and 
Type  B  S3^stems. 

For  the  8  Weibull  components  in  parallel-series  (Figure  4.3)  in  Case 
8,  no  change  was  made  for  the  mission  time  for  each  of  the  components 
which  is  10  hours.  The  ratios  between  the  largest  and  the  smallest  failure 
rate  were  chosen  to  be  8  for  both  systems.  The  shape  parameter,  j?;,  was 
set  at  2  for  all  subcases. 

Case  9  is  a  mixture  of  exponential  and  Weibull  components  in 
parallel-series  (Figure  4.3).  The  mission  time  for  all  components  is  10 
hours.  The  values  of  the  scale  parameters,  /l„  for  exponential  compo¬ 
nents  were  chosen  to  be  equal  to  those  for  Weibull  components.  The 
ratios  between  the  largest  and  the  smallest  failure  rate  were  chosen  to  be 
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4  for  both  component  types  and  both  systems.  The  shape  parameter,  pi, 
was  set  at  2  for  all  subcases. 

Case  10  is  a  mixture  of  exponential  and  Weihull  components  in  a 
more  general  series-parallel  configuration  (Figure  4.4).  The  mission  time 
for  each  of  the  components  was  chosen  randomly.  Similarly,  the  scale 
parameter  for  each  of  the  components  and  shape  parameter  for  each  of 
the  Weibull  components  were  chosen  randomly.  This  was  done  for  both 
systems. 

Each  simulation  run  of  1000  replications  results  in  an  output  file 
OUT.DAT.  The  raw  output  from  all  the  RETP  runs  are  summarized  in 
tabular  form  and  placed  in  Appendix  E.  Each  Table  from  Table  lA  to 
Table  3B  corresponds  to  a  specific  run  case  and  system  type  combination 
using  various  degrees  of  freedom.  Each  Table  from  Table  4A  to  Table 
lOB  corresponds  to  a  specific  run  case  and  system  type  combination  using 
only  one  specific  degree  of  freedom. 
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V.  ANALYSIS  OF  SIMULATION  RESULTS 


After  simulation  runs  were  completed  for  selected  cases,  the  results 
were  analyzed  for  comparisons  with  results  that  had  been  obtained  by 
Yee  [Ref,  7]  . 

For  each  case,  simulations  were  run  with  different  degrees  of  freedom 
for  the  chi-square  percentile  point.  This  was  done  to  determine  if  a  for¬ 
mula  for  the  degrees  of  freedom  could  be  developed  that  would  5deld  a 
more  accurate  lower  confidence  limit  procedure. 

TEST  PLAN :  TESTING  COMPONENTS  UNTIL  fi  FAIL  (RETP) 

When  all  components  of  the  system  have  exponentially  distributed 
failure  times,  the  lower  confidence  limit  procedures  in  this  thesis  are 
nearly  the  same  as  those  developed  by  Yee.  Consequently  analysis  of  the 
simulation  results  will  be  discussed  primarily  for  cases  when  some  com¬ 
ponents  of  the  system  have  failure  times  that  have  a  Weibull  distribution. 

A  comparison  of  the  four  values  of  RSLOW,  as  in  Table  lA  for  each 
of  5  sampling  plans  (denoted  by  S/N),  reveals  that  the  lower  confidence 
limit  procedure  with  degrees  of  freedom  equal  to  2NFC  -  NCOMP  is  the 
most  accurate  lower  confidence  limit  procedure.  In  S/N  2,  for  example, 
the  RSLOW  value  of  0.9298  (using  2NFC  -  NCOMP  degrees  of  freedom) 
is  the  largest  such  value  below  the  RS  value  of  0.9305.  The  value  of 
RSLOW  above  RS  are  optimistic  and  not  as  desirable  as  values  of 
RSLOW  which  are  equi-distant  below  RS. 

Table  2A  and  Yee's  Table  2A  provide  a  comparison  of  the  simulation 
results  for  eight  Weibull  components  in  series  using  the  same  key  pa¬ 
rameters  and  the  same  mission  times  for  each  component  as  used  by  Yee. 
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When  testing  15  components  until  all  fail,  Yee's  procedure  and  the 
procedure  in  this  thesis  have  nearly  equal  accuracy.  But  when  testing  15 
items  until  7  fail,  the  procedure  written  in  this  thesis  gives  a  more  accu¬ 
rate  value  of  RSLOW  than  the  one  from  Yee's  procedure.  In  this  case, 
Yee's  procedure  yields  values  of  RSLOW  that  are  all  above  RS.  In  both 
Yee's  procedure  and  the  one  used  in  this  thesis,  the  degrees  of  freedom 
of  2(NFC  -  NCOMP)  gave  the  most  accurate  results  among  the  four 
choices  of  degrees  of  freedom.  It  can  also  be  seen  that  for  more  truncated 
testing,  such  as  testing  15  items  until  3  fail,  the  value  of  RSLOW  tends 
to  be  higher  and  they  are  all  slightly  above  RS. 

Table  3 A  displays  the  results  of  Case  3  for  a  type  A  system.  The  shape 
parameter,  j5,  for  the  failure  times  with  Weibull  distribution  was  set  equal 
to  2.  The  value  of  RSLOW  resulting  from  the  procedure  developed  in 
this  thesis  is  more  accurate  than  that  for  Yee's  procedure  for  various  de¬ 
grees  of  freedom.  Inspection  of  Table  3A  reveals  that  the  procedure 
corresponding  to  degrees  of  freedom  2(NFC  f  NCOMP)  is  reasonably 
accurate  for  all  5  simulation  cases  and  for  both  90%  and  80%  confidence 
levels. 

Since  the  accuracy  of  the  procedure  depends  on  the  extent  of  the 
truncation  in  the  testing,  a  method  was  developed  for  choosing  the  de¬ 
grees  of  freedom  that  includes  a  term  with  the  ratio  NC/NF,  where  NC 
is  the  number  of  components  placed  on  test  and  NF  is  the  number  of 
failed  components.  This  makes  the  procedure  dependent  on  the  amount 
of  truncation  in  the  testing.  Using  information  from  all  the  simulations 
and  the  degrees  of  freedom  that  gave  the  most  accurate  values  of  RSLOW 
for  all  cases  simulated,  a  formula  for  the  degrees  of  freedom,  DF,  was 
developed  using  least  squares  methods. 
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The  equation  is  DF  =  2NFC  +  1.1(NC/NF)  (see  Appendix  B  for  the 
derivation  of  this  equation).  Applying  the  equation  to  obtain  RSLOW 
and  observing  the  results  for  all  cases  simulated,  resulted  in  a  small 
modification  to  the  above  equation  to  yield  the  final  equation  DF  = 
2NFC  +  0.5(NC/NF).  This  final  equation  was  used  for  the  remaining 
cases  that  were  simulated  in  RETP  when  any  components  had  failure 
times  with  Weibull  distribution. 
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Table  lA  :  8  Exp  in  Series,  RS  =  0.9305  (Hi) 


min  A  =  C 

1.0002  f/hr,  max  A  = 

0.0016  f/hr,  UT  =  10  hrs. 

S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

oc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.9262 

0.9630 

10  failed 

(160) 

0.2 

0.9256 

0.9220 

NFC  =80 

2(NFC  + 

0.1 

0.9196 

0.9970 

NCOMP) 

(176) 

0.2 

0.9188 

0.9910 

2NFC- 

0.1 

0.9296 

0.9209 

NCOMP 

(152) 

0.2 

0.9291 

0.8430 

2(NFC- 

0.1 

0.9329 

0,8400 

NCOMP) 

(144) 

0.2 

0.9325 

0.7300 

2 

Test  15  until 

2NFC 

0.1 

0.9277 

0.9550 

15  failed 

(240) 

0.2 

0.9273 

0.9080 

NFC=  120 

2(NFC+ 

0.1 

0.9233 

0.9899 

NCO.MP) 

(256) 

0.2 

0.9228 

0.9750 

2NFC- 

0.1 

0.9298 

0.9159 

NCOMP 

(232) 

0.2 

0.9296 

0.8329 

2(NFC- 

0.1 

0.9321 

0.8440 

NCOMP) 

(224) 

0.2 

0.9318 

0.7470 

3 

Test  15  until 

2NFC 

0.1 

0.9268 

0.9550 

11  failed 

(176) 

0.2 

0.9262 

0.9159 

NFC  =88 

2(NFC+ 

0.1 

0.9208 

0.9960 

NCOMP) 

(192) 

0.2 

0.9200 

0.9880 

2NFC- 

0.1 

0.9298 

0.9159 

NCOMP 

(168) 

0.2 

0.9293 

0.8430 

2(NFC- 

0.1 

0.9328 

0.8430 

NCOMP) 

(160) 

0.2 

0.9324 

0.7350 
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Table  lA  :  8  Exp  in  Series,  RS  =  0.9305 
min  ).  =  0.0002  f/hr,  max  X 


(Hi)  (Com...) 


Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

cc 

RSLOW 

LEVEL 

Test  15  until 

2NFC 

0.1 

0.9241 

0.9700 

7  failed 

(112) 

0.2 

0.9229 

0.9310 

NFC =56 

2(NFC  + 

0.1 

0.9145 

0.9980 

NCOMP) 

(128) 

0.2 

0.9130 

0.9940 

2NFC- 

0.1 

0.9289 

0.9190 

NCOMP 

(104) 

0.2 

0.9280 

0.8530 

2(NFC- 

0.1 

0.9338 

0.8350 

NCOMP) 

(96) 

0.2 

0.9331 

0.7200 

Test  15  until 

2NFC 

0.1 

0.9145 

0.9859 

3  failed 

(48) 

0.2 

0.9129 

0.9750 

NFC =24 

2(NFC+ 

0.1 

0.8907 

1.0000 

NCOMP) 

(64) 

0.2 

0.8877 

1.0000 

2NFC- 

0.1 

0.9268 

0.9439 

NCOMP 

(40) 

0.2 

i 

0.9260 

0.8600 

2{NFC- 

0.1 

0.9394 

0.7530 

NCOMP) 

(32) 

0.2 

0.9394 

0.6339 
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Yee's  resuits 


Table  lA  :  S  Exp  in  Series,  RS  =  0.931  (Hi) 

_  min  X  -  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

.Measures  of  accuracy 

RSLOW 

LEVEL 

I 

Test  5  until 

5  Tailed 

NFC =40 

2NFC 

(80) 

0.1 

0.919 

0.982 

0.2 

0.919 

0.960 

2{NFC+ 

NCOMP) 

(r6) 

O.i 

0.906 

1.000 

0.2 

0.905 

0.999 

2NFC- 

NCOxMP 

(72) 

0.1 

0.927 

0.949 

0.2 

0.927 

0.880 

2(NFC. 

NCOMP) 

(64) 

0.1 

0.934 

0.821 

0.2 

0.934 

0.702 

2 

Test  15  until 

15  failed 

NFC  =120 

2NFC 

(240) 

0.1 

0.928 

0.955 

0.2 

0.927 

0.908 

2(NFC+ 

NCOMP) 

(256) 

0.1 

0.923 

0.950 

0.2 

0.923 

0.975 

2NFC. 

NCOMP 

(232) 

0.1 

0.930 

0.916 

0.2 

0.930 

0.833 

2(NFC- 

NCO.MP) 

(224) 

0.1 

0.53'’ 

0.844 

0.2 

0.932 

0.747 

3 

Test  15  until 

11  failed 

NFC  =88 

2NFC 

(176) 

0.1 

0.927 

0.955 

0.2 

0.926 

0.916 

2(NFC  + 
NCOMP) 
(192) 

0.1 

0.921 

0.996 

0.2 

0.920 

0.988 

2NFC- 

NCO.MP 

(168) 

0.1 

0.930 

0.916 

0.2 

0.92« 

0.843 

2(NFC- 

NCOMP) 

(160) 

0.1 

0.933 

0.843 

0.2 

0.932 

0.735 
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Yee's  results 


Table  lA  :  8  Exp  in  Series,  RS  =  0.931  (Hi)  (Cont...) 

_  min  k  —  0.0002  f/hr,  max  A  =  0.0016  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Mfasures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC =56 

2NFC 

0.1 

0.924 

0.970 

(112) 

0.2 

0.923 

0.931 

2{NFC+ 

NCOMP) 

(;28) 

0.1 

0.915 

0.998 

0.2 

0.913 

0.994 

2NFC- 

NCOMP 

(104) 

0.1 

0.929 

0.919 

0.2 

0.928 

0.853 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.934 

0.835 

0.2 

0.933 

0.720 

5 

Test  15  until 

3  failed 

NFC =24 

2NFC 

(48) 

0.1 

0.915 

0.986 

0.2 

0.912 

0.975 

2(NFC+ 

NCOMP) 

(64) 

0.1 

0.891 

1.000 

0.2 

0.888 

1.000 

2NFC- 

NCOMP 

(40) 

0.1 

0.927 

0.944 

0.2 

0.926 

0.860 

2(NFC- 

NCO.MP) 

(32) 

0.1 

0.939 

0.753 

0,2 

0.939 

0.634 
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Table  2A  :  8  Wei  in  Series,  RS  =  0.9798  (Hi) 

_ min  ).  =  0.001  f/hr,  max  A  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

cc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.9602 

0.9859 

10  failed 

(160) 

0.2 

0.9506 

0.9820 

NFC  =  80 

2(NFC+ 

0.1 

0.9565 

0.9890 

NCOMP) 

(176) 

0.2 

0.9461 

0.9870 

2NFC. 

0.1 

0.9619 

0.9820 

NCOMP 

(152) 

0.2 

0.9529 

0.9809 

2(NFC- 

0.1 

0.9638 

0.9800 

NCOMP) 

(144) 

0.2 

0.9553 

0.9770 

2 

Test  15  until 

2NFC 

0.1 

0.9706 

0.9770 

15  failed 

(240) 

0.2 

0.9649 

0.9719 

NFC  =120 

2(NFC+ 

0.1 

0.9687 

0.9820 

] 

NCOMP) 

(256) 

0.2 

0.9627 

0.9790 

2NFC- 

0.1 

0.9715 

0.9730 

NCOMP 

(232) 

0.2 

0.9659 

0.9640 

2(NFC- 

0.1 

0.9724 

0.9650 

NCOMP) 

(224) 

0.2 

0.9671 

0.9590 

3 

Test  15  until 

2NFC 

0.1 

0.9704 

0.9719 

11  failed 

(176) 

0.2 

0.9637 

0.9650 

NFC  =88 

2(NFC+ 

0.1 

0.9679 

0.9800 

NCOMP) 

(192) 

0.2 

0.9606 

0.9760 

2NFC-  1 

0.1 

0.9716 

0.9660 

j 

NCOMP 

(168) 

0.2 

0.9652 

0.9590 

2(NFC- 

0.1 

1 

0.9729 

0.9590 

NCOMP) 

(160) 

0.2 

0.9668 

0.9510 
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Table  2A  :  8  Wei  in  Series,  RS  =  0.9798  (Hi)  (Cent...) 


min  A  ==  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

.Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC  =56 

2NFC 

(112) 

0.1 

0.9764 

0.9400 

0.2 

0.9668 

0.9280 

2(NFC  + 
NCOMP) 
(128) 

0.1 

0.9733 

0.9579 

0.2 

0.9624 

0.9510 

2NFC- 

NCO.MP 

(104) 

0.1 

0.9779 

0.9240 

0.2 

0.9691 

0.9110 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.9795 

0.9080 

0.2 

0.9713 

0.8900 

5 

Test  15  until 

3  failed 

NFC  =24 

2NFC 

(48) 

0.1 

0.9889 

0.8000 

0.2 

0.9814 

0.7770 

2(NFC  + 
NCO.MP) 
(64) 

0.1 

0.9857 

0.8469 

0.2 

0.9758 

0.8360 

2NFC- 

NCOMP 

(40) 

0.1 

0.9906 

0.7510 

0.2 

0.9843 

0.7280 

2(NFC. 

NCOMP) 

(32) 

0.1 

0.9922 

0.6890 

0.2 

0.9872 

0.6530 
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Yee's  results 

Table  2A  :  8  Wei  in  Series,  RS  =  0.980  (Hi) 

_ min  X  =  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

CL 

RSLOW 

LEVEL 

1 

Test  5  until 

2NFC 

0.1 

0.947 

0.992 

5  failed 

(80) 

0.2 

0.930 

0.989 

NFC  =  40 

2(NFC+ 

0.1 

0.937 

0.994 

NCOMP) 

(96) 

0.2 

0.918 

0.993 

2NFC- 

0.1 

0.951 

0.989 

NCOMP 

(72) 

0.2 

0.937 

0.986 

2(NFC- 

0.1 

0.956 

0.985 

NCOMP) 

(64) 

0.2 

0.943 

0.981 

2 

Test  15  until 

2NFC 

0.1 

0.978 

0.918 

15  failed 

(240) 

0.2 

0.974 

0.913 

NFC  =120 

2(NFC+ 

0.1 

0.977 

0.931 

NCO.MP) 

(256) 

0.2 

0.972 

0.924 

2NFC- 

0.1 

0.979 

0.914 

NCOMP 

(232) 

0.2 

0.975 

0.901 

2(NFC- 

0.1 

0.980 

0.904 

NCO.MP) 

(224) 

0.2 

0.975 

0.889 

3 

Test  15  until 

2NFC 

0.1 

0.982 

0.876 

11  failed 

(176) 

0.2 

0.977 

0.860 

NFC  =88 

2(NFC+ 

0.1 

0.980 

0.894 

NCOMP) 

(192) 

0.2 

0.975 

0.882 

2NFC- 

0.1 

0.983 

0.861 

NCOMP 

(168) 

0.2 

0.978 

0.839 

2(NFC- 

0.1 

0.983 

0.840 

NCOMP) 

(160) 

0.2 

0.979 

0.819 
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Yee's  results 

Table  2A  :  8  Wei  in  Series,  RS  =  0.980  (Hi)  (Cent...) 


min  X  =  0.001  f/hr,  max  ).  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC  =  56 

2NFC 

(112) 

0.1 

0.987 

0.800 

1  0.2 

0.981 

0.779 

2(NFC  + 
NCOMP) 
(128) 

0.1 

0.985 

0.839 

0.2 

0.978 

0.824 

2NFC- 

NCOMP 

(104) 

0.1 

0.988 

0.776 

0.2 

0.982 

0.753 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.989 

0.746 

0.2 

0.983 

0.732 

5 

1 

1 

! 

1 

Test  15  until 

3  failed 

NFC  =24 

1 

2NFC 

(48) 

0.1 

0.994 

0.621 

0.2 

0.991 

0.584 

2(NFC  + 
NCOMP) 
(64) 

0.1 

0.993 

0.705 

0.2 

0.988 

0.685 

2NFC- 

NCOMP 

(40) 

0.1 

0.995 

0.548 

0.2 

0.992 

0.514 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.996 

0.408 

0.2 

0.993 

i 

0.417 
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Table  3A  :  4  Exp  and  4  Wei  (mixed)  in  Series,  RS  =  0.9801  (Hi) 
_ min  ).  =  0.0002  f/hr,  max  X  —  0.0008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

oc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

o.i 

0.9796 

0.9190 

10  failed 

(160) 

0.2 

0.9791 

0.8720 

NFC  =80 

2(NFC+ 

0.1 

0.9777 

0.9740 

NCOMP) 

(176) 

0.2 

0.9771 

0.9489 

2NFC- 

0.1 

0.9805 

0.8720 

\ 

NCO.MP 

(152) 

0.2 

0.9801 

0.7980 

2(NFC- 

0.1 

0.9815 

0.7950 

NCOMP) 

(144) 

0.2 

0.9811 

0.7000 

2 

Test  15  until 

2NFC 

0.8950 

15  failed 

(240) 

0.9799 

0.8250 

NFC=  120 

2(NFC  + 

0.1 

0.9789 

0.9510 

NCOMP) 

(256) 

0.2 

0.9786 

0.9209 

2NFC- 

0.1 

0.9808 

0.8370 

NCOMP 

(232) 

0.2 

0.9805 

0.7530 

2(NFC- 

0.1 

0.9814 

0.7670 

NCOMP) 

(224) 

0.2 

0.9812 

0.6670 

3 

Test  15  until 

2NFC 

0.1 

0.9802 

0.8920 

11  failed 

(176) 

0.2 

0.9797 

0.8400 

NFC  =88 

2(NFC+ 

0.1 

0.9785 

0.9590 

NCOMP) 

(192) 

0.2 

0.9779 

0.9320 

2NFC- 

0.1 

0.9810 

0.8410 

NCOMP 

(168) 

0.2 

0.9806 

0.7550 

2(NFC- 

0.1 

0.9819 

0.7580 

NCOMP) 

(160) 

0.2 

0.9815 

0.6590 
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Table  3A  :  4  Exp  and  4  Wei  in  Series,  RS  =  0.9801  (Hi)  (Cent...) 


min  X  =  0.0002  f/hr,  max  A  =  0.0008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC  =  56 

2NFC 

(112) 

0.1 

0.9803 

0.8950 

0.2 

0.9797 

0.8260 

2(NFC  + 
NCOMP) 
(128) 

0.1 

0.9777 

0.9690 

0.2 

0.9770 

0.9439 

2NFC. 

NCO.MP 

(104) 

0.1 

0.9815 

0.8060 

0.2 

0.9811 

0.7300 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.9828 

0.7090 

0.2 

0.9825 

0.6000 

5 

Test  15  until 

3  faded 

NFC  =24 

2NFC 

(48) 

0.1 

0.9815 

0.8730 

0.2 

0.9800 

0.8020 

2(NFC  + 
NCOMP) 
(64) 

0.1 

0.9757 

0.9740 

0.2 

0.9739 

0.9560 

2NFC- 

NCOMP 

(40) 

0.1 

0.9839 

0.7250 

0.2 

0.9831 

0.6260 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.9868 

0.5150 

0.2 

0.9862 

0.4140 

Yee's  results 

Table  3A  :  4  Exp  and  4  Wei  i^mixed)  in  Series,  RS  =  0.980  (Hi) 
_ min  ).  =  0.002  f/hr,  max  A  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

cc 

RSLOW 

LEVEL 

1 

Test  5  until 

2NFC 

0.1 

0.979 

0.942 

5  failed 

(80) 

0.2 

0.978 

0.905 

NFC =40 

2(NFC+ 

0.1 

0.975 

0.987 

NCOMP) 

(96) 

0.2 

0.974 

0.976 

2NFC- 

0.1 

0.981 

0.881 

NCOiMP 

(72) 

0.2 

0.980 

0.805 

2(NFC- 

0.1 

0.983 

0.771 

NCOMP) 

(64) 

0.2 

0.982 

0.684 

2 

Test  15  until 

2NFC 

0.1 

0.981 

0.863 

15  failed 

(240) 

0.2 

0.980 

0.800 

NFC  =120 

2(NFC+ 

0.1 

0.979 

0.941 

NCOMP) 

(256) 

0.2 

0.979 

0.898 

2NFC- 

0.1 

0.981 

0.881 

NCOMP 

(232) 

0.2 

0.980 

0.805 

2(NFC. 

0.1 

0.982 

0.725 

NCOMP) 

(224) 

0.2 

0.981 

0.631 

3 

Test  15  until 

2NFC 

0.1 

0.981 

0.864 

11  failed 

(176) 

0.2 

0.980 

0.801 

1 

NFC  =  88 

2(NFC  + 

0.1 

0.979 

0.951 

i 

NCOMP) 

(192) 

0.2 

0.978 

0.907 

2NFC- 

0.1 

0.982 

0.802 

j 

NCO.MP 

(168) 

0.2 

0.981 

0.698 

2(NFC- 

0.1 

0.982 

0.702 

NCOMP) 

(160) 

0.2 

0.982 

0.591 

39 


Yee's  results 


Table  3A  :  4  Exp  and  4  Wei  in  Series,  RS  =  0.980  (Hi)  (Cont...) 
_ min  ).  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC =56 

2NFC 

(112) 

0.1 

0.981 

0.865 

0.2 

0.980 

0.787 

2(NFC  + 
NCOMP) 
(128) 

0.1 

0.978 

0.952 

0.2 

0.978 

0.920 

2NFC- 

NCO.MP 

(104) 

0.1 

0.982 

0.769 

0.2 

0.982 

0.676 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.983 

0.644 

0.2 

0.983 

0.523 

5 

Test  15  until 

3  failed 

NFC =24 

2NFC 

(48) 

0.1 

0.982 

0.843 

0.2 

0.981 

0.762 

2(NFC+ 

NCOMP) 

(64) 

0.1 

0.976 

0.970 

0.2 

0.975 

0.941 

2NFC- 

NCOMP 

(40) 

0.1 

0.984 

0.684 

0.2 

0.984 

0.580 

2(NFC- 

NCO.MP) 

(32) 

0.1 

0.987 

•0.459 

0.2 

0.987 

0.386 

Table  4A  :  8  Exponential  in  Parallel,  RS  =  0.9345  (Hi) 

min  X  =  0.1000  f/hr,  max  X  =  0.1525  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

160 

0.1 

0.9306 

0.9290 

0.2 

0.9320 

0.8180 

2 

Test  15  until 

15  failed 

NFC  =120 

240 

0.1 

0.9325 

0.9190 

0.2 

0.9324 

0.8210 

3 

Test  15  until 

11  failed 

NFC  =88 

176 

0.1 

0.9309 

0.9240 

0.2 

0.9312 

0.8329 

4 

Test  15  until 

7  failed 

NFC =56 

112 

0.1 

0.9307 

0.9209 

0.2 

0.9308 

0.8310 

5 

Test  15  until 

3  failed 

NFC  =24 

48 

0.1 

0.9280 

0.9220 

0.2 

0.9286 

0.8370 
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Table  5A  :  8  Wei  in  Parallel  ,  RS  =  0.9265  (Hi) 

min  X  =  0.100  f/hr,  max  X  -  0.128  f/hr,  UT  =  10  hrs. 


Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.9154 

0.9540 

0.2 

0.9157 

0.8880 

2 

Test  15  until 

15  failed 

NFC=  120 

241 

0.1 

0.9179 

0.9510 

0.2 

0.9170 

0.8770 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9248 

0.9180 

0.2 

0.9250 

0.8130 

4 

Test  15  until 

7  failed 

NFC  =  56 

113 

0.1 

0.9323 

0.8630 

0.2 

0.9318 

0.7589 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.9152 

0.9270 

0.2 

0.8980 

0.8850 
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Table  6A  :  4  EXP  and  4  Wei  (mixed)  in  Parallel,  RS  =  0.9408  (Hi) 
min  X  =  0.100  f/hr,  max  X  =  0.130  f/hr,  UT  =  10  hrs. 


Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

oc 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.9342 

0.9420 

0.2 

0.9366 

0.8500 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.9362 

0.9330 

0.2 

0.9363 

0.85)0 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9389 

0.9240 

g 

0.9378 

0.8310 

4 

Test  15  until 

7  failed 

NFC  =  56 

113 

m 

0.9416 

0.8920 

0.2 

0.9427 

0.7819 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.9383 

0.9080 

0.2 

0.9325 

0.8410 

Table  7A  :  8  Exp  in  Series-  Parallel,  RS  =  0.9249  (Hi) 

min  X  =  0.0003  f/hr,  max  X  =  0.0024  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

160 

0.1 

0.9241 

0.9159 

0.2 

0.9226 

0.8620 

2 

Test  15  until 

15  failed 

NFC  =120 

240 

0.1 

0.9253 

0.8960 

0.2 

0.9236 

0.8360 

3 

Test  15  until 

11  failed 

NFC  =88 

176 

0.1 

0.9252 

0.8950 

0.2 

0.9226 

0.8390 

4 

Test  15  until 

7  failed 

NFC =56 

112 

0.1 

0.9222 

0.9260 

0.2 

0.920^  1  “  >>  •  ,9 

5 

Test  15  until 

3  failed 

NFC =24 

48 

0.1 

0.2 

» 

0.9123  1  O.9260 

_ 1 _ 
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Table  8A  :  8  Wei  in  Series-  Parallel,  RS  =  0.9328  (Hi) 

min  X  =  0.002  f/hr,  max  X  =  0.016  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  n-  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.9106 

0.9579 

0.2 

0.8939 

0.9510 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.9247 

0.9330 

0.2 

0.9108 

0.9230 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9266 

0.9220 

0.2 

0.9149 

0.9119 

4 

Test  15  until 

7  failed 

NFC  =  56 

113 

0.1 

0.9414 

0.8609 

0.2 

WM 

0.8419 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.7079 

0.2 

0.9536 

0.6740 
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Table  9A  :  4  EXP  and  4  Wei  in  Series-  Parallel,  RS  =  0.9276  (Hi) 
min  X  =  0.005  f/hr,  max  X  -  0.020  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =  80 

161 

0.1 

0.9096 

0.9400 

0.2 

0.8904 

0.9349 

■  2 

Test  15  until 

15  failed 

NFC=  120 

241 

0.1 

0.9238 

0.9170 

0.2 

0.9079 

0.9040 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9266 

0.9050 

0.2 

0.9109 

0.8850 

4 

Test  15  until 

7  failed 

NFC =56 

113 

0.1 

0.9393 

0.8550 

0.2 

0.9203 

0.8390 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.9637 

0.7430 

0.2 

0.9451 

0.7200 
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Table  lOA  :10  EXP  and  5  Wei  in  Series-  Parallel,  RS  =  0.9472  (Hi) 
Exp  :  min  X  =  0.025  f/hr,  max  X  =  0.075  f/hr. 

Wei  :  min  X  =  0.055  f/hr,  max  X  =  0.095  f/hr. 

UT  =  10  hrs. 


Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

I 

Test  10  until 

10  failed 

NFC  =150 

301 

0.1 

0.9383 

0.9710 

0.2 

0.9369 

0.9550 

2 

Test  20  until 

20  failed 

NFC=  300 

601 

0.1 

0.9449 

0.9370 

0.2 

0.9432 

0.8950 

3 

Test  20  until 

15  failed 

NFC =225 

451 

0.1 

0.9444 

0.9349 

0.2 

0.9430 

0.8979 

4 

Test  20  until 

10  failed 

NFC=  150 

301 

0.1 

0.9435 

0.9290 

0.2 

0.9415 

0.9000 

5 

_ ... 

Test  20  until 

5  failed 

NFC  =75 

152 

0.1 

0.9344 

0.9809 

0.2 

0.9310 

0.9520 
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VI.  CONCLUSIONS 


The  lower  confidence  limit  procedures  developed  in  this  thesis  are 
modifications  and  extensions  of  the  lower  confidence  limit  procedures 
written  by  Yee.  Procedures  were  developed  and  evaluated  in  this  thesis  for 
more  complex  system  structures  than  Yee  examined.  In  addition  a  dif¬ 
ferent  method  for  estimating  the  parameters,  in  the  Weibull  distribution 
was  used  here  rather  than  the  maximum  likelihood  procedure  used  by 
Yee. 

The  evaluations  show  that  the  approximate  lower  confidence  limit 
procedures  are  reasonably  accurate  for  all  system  structures  examined  if 
the  degrees  of  freedom  are  chosen  judiciously.  Although  the  equations 
given  here  for  choosing  an  appropriate  value  for  the  degrees  of  freedom 
are  for  the  system  simulated  in  this  thesis,  it  would  be  prude  it  to  run 
simulations  for  complex  systems  that  differ  substantially  from  those  ex¬ 
amined  in  this  thesis  in  order  to  determine  an  appropriate  number  for  the 
degrees  of  freedom  in  the  lower  confidence  limit  equation. 

The  degrees  of  freedom  equation  derived  in  Appendix  B,  DF  = 
2NFC  +  0.5(NC/NF),  yielded  lower  confidence  limit  procedures  with 
reasonable  accuracy  for  test  plans  that  were  truncated  on  the  number  of 
failures. 
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VII.  APPLICATION  EXAMPLES 


Based  on  the  procedures  evaluated  by  the  RETP  runs,  three  config¬ 
urations  of  systems,  a  specific  test  plan  and  failure  time  data  sets  were 
constructed  to  illustrated  the  use  of  the  procedures  in  providing  a  lower 
100(1  —  a)  %  confidence  limit  for  system  reliability.  Actual  results  of  one 
computer  run  for  each  case  simulated  are  given  below. 


Case  1  ;  8  Exponential  components  in  Series 

TEST  PLAN  1  -  Test  15  tmtil  7  fails  for  each  component 

I.  Raw  Data 


Comp  Ordered  Failure  Times 


1 

T(l) 

T(2) 

TC3) 

T(4) 

TC5) 

T(6) 

T(7) 

1 

695.0 

1241.  3 

1365. 7 

2628. 8 

3304. 9 

3946. 3 

4014.  r: 

2 

142.5 

212.7 

230.8 

315.2 

401.4 

1071.9 

1222. 0 

3 

111.9 

394.7 

422.6 

506.  5 

519.5 

558.5 

582.6 

4 

325.5 

356.2 

441.7 

837.  8 

844.2 

873.7 

894.4 

5 

62.6 

110.0 

124.3 

126.  8 

325.8 

384.0 

502.  8 

6 

20.0 

85.6 

107.0 

108.  0 

161.5 

190. 1 

201.6 

7 

34. 1 

77.9 

100.4 

142.  1 

156.0 

180.9 

193.5 

8 

6.1 

34.6 

48.  1 

65.  3 

95.2 

136.6 

160.4 

II.  Data  Summary 


CfjmpCi)  UT(i) 

NC(1) 

NF(i) 

TTCi) 

ELM(i) 

ER(i) 

ER(i)*TT(i) 

1 

10.0 

15 

7 

49310.2 

0.00012 

0. 03711 

1829.9 

2 

10.0 

15 

7 

13372.8 

0. 00045 

0.  13683 

1829.  9 

3 

10.0 

15 

7 

7757.2 

0.00077 

0.23589 

1829. 9 

4 

10.0 

15 

7 

11728.  8 

0. 00051 

0. 15601 

1829. 9 

5 

10.0 

15 

7 

5658.4 

0.00106 

0. 32339 

1829.9 

6 

10.0 

15 

7 

2487. 0 

0.  00241 

0.  73578 

1829. 9 

7 

10.0 

15 

7 

2433. 2 

0.  00247 

0. 75205 

1829. 9 

8 

10.  0 

15 

7 

1829. 9 

0. 00328 

1.00000 

1829. 9 

III.  Estimation  Procedure  for  RSLOW 


Parameter  df  Value 
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RS  0,530r'. 

Ar-r;  o.  i 

NFC  56 

CHISQD  104  122.86  (from  the  table  Chi-square  distribution  ) 

LMU  0. 00419 


RSLOW  0. 92894 


CASE  4  :  8  Exponential  components  in  parallel 
- -  TEST  PLAN  1  -  Test  15  mtil  7  fails  for  each  component 

I.  Raw  Data 


Comp  Ordered  Failure  Times 


i 

T(l) 

TC2) 

T(3) 

T(4) 

T(5) 

T(6) 

T(7) 

1 

1. 3899 

2. 4825 

2.7315 

5.2576 

6.6099 

7.8927 

8. 0285 

2 

0.5303 

0.7915 

0. 8586 

1. 1727 

1.4936 

3.9886 

4.5471 

3 

0.5839 

2.  0591 

2. 2047 

2. 6425 

2.7106 

2.9138 

3.0398 

4 

2. 1258 

2.  3264 

2. 8844 

5.4715 

5.5130 

5.7055 

5.8410 

5 

0.4813 

0. 8465 

0.9563 

0.9751 

2.5065 

2. 9542 

3.8674 

6 

0. 1748 

0.7471 

0. 9340 

0. 9424 

1.4090 

1.6593 

1.7597 

7 

0.3296 

0.  7519 

0.9697 

1.3716 

1.5065 

1. 7469 

1.8685 

8 

0. 0638 

0. 3635 

0.5050 

0, 6849 

0.9989 

1.  4334 

1. 6832 

II.  Data  Summary 


Comp(i)  UT(i) 

NC(i) 

NF(i) 

TT(i) 

ELM(i) 

ER(i) 

ER(i)-Tr(i) 

1 

10.  0 

15 

7 

98. 6204 

0.  06084 

0.  19467 

19. 1986 

2 

10.0 

15 

7 

49. 7592 

0.  12058 

0.  38583 

19. 1986 

3 

10.0 

15 

7 

40.4724 

0.  14825 

0.47436 

19. 1986 

4 

10.0 

15 

7 

76.5961 

0.07833 

0.  25065 

19. 1986 

5 

10.0 

15 

7 

43.5263 

0.  13785 

0.  44108 

19. 1986 

6 

10.0 

15 

7 

21.  7044 

0.27644 

0.  88455 

19.  1986 

7 

10.0 

15 

7 

23.4926 

0.  25540 

0.81722 

19.  1986 

8 

10.0 

15 

7 

19.1986 

0.31252 

1.00000 

19. 1986 

III.  Estimation  Procedure  for  RSLOW 


irameter 

df 

Value 

RS 

0.  93459 

ALPHA 

0. 1 

NFC 

56 

CHISQD 

112 

131.56  (from  the  table  Chi-square  distribution  ) 

LMU 

0.42828 
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RSLOW 


0.93069 


Case  7  :  8  Exponential  components  in  Series-Parallel. 

-  (4  components  in  parallel  connected  to  4  components  in  series  ) 

TEST  PLAN  1  -  Test  15  until  7  fails  for  each  component 

I.  Raw  Data 


Comp  Ordered  Failure  Times 


i 

T(l) 

T(2) 

T(3) 

T(4) 

T(5) 

TC6) 

T(7) 

1 

463. 309 

827.505 

910.497 

1752.519 

2203. 295 

2630. 893 

2676. 161 

2 

95.019 

141.  817 

153. 837 

210. 102 

267. 609 

714.  622 

814.687 

3 

74. 604 

263. 110 

281.  713 

337. 648 

346.  349 

372.313 

388.416 

4 

217.008 

237.491 

294.  454 

558.552 

562.  782 

582.436 

596.273 

5 

41.  716 

73. 360 

82. 882 

34.511 

217.230 

256. 028 

335. 172 

6 

13. 354 

57. 073 

71.346 

71.990 

107. 635 

126. 749 

134.425 

7 

22. 755 

51  917 

66. 954 

94. 703 

104. 022 

120.616 

129.016 

8 

4.  053 

23. 095 

32.  090 

43.517 

63.475 

91.  080 

106. 955 

II.  Data  Summary 


Comp(i)  UT(i) 

NCCi) 

NF(i) 

TT(i) 

ELM(i) 

ERCi) 

FRCi)*TT(i) 

1 

10.0 

15 

7 

32873.465 

0.  00018 

0. 03711 

1219.9089 

2 

10.0 

IS 

/ 

8915. 184 

0. 00067 

0.  13683 

1219.9089 

3 

10,  U 

15 

7 

5171.477 

0. 00116 

0.  23589 

1219.  9089 

4 

10.0 

15 

7 

7819. 181 

0.00077 

0.  15601 

1219.  9089 

5 

10.0 

15 

7 

3772.275 

0.  00159 

0.  32339 

1219.  9089 

6 

10,  0 

15 

7 

1657.971 

0.  00362 

0.  73578 

1219.9089 

7 

10.0 

15 

7 

1622. 109 

0.  00370 

0.  75205 

1219.9089 

8 

10.  0 

15 

7 

1219.909 

0.  00492 

1.  00000 

1219. 9089 

III.  Estimation  Procedure  for  RSLOW 


Parameter 

df 

Value 

RS 

0.  92496 

ALPHA 

0.  1 

NTC 

56 

CHISQD 

112 

131.56  (from 

LMU 

0.00674 

RSLOW 

0.  92223 

the  table  Chi-square  distribution  ) 


Case  2  :  8  Weibull  components  in  Series 

- - —  TEST  PLAN  1  -  Test  15  xmtil  7  fails  for  each  component 
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1.  Raw  Data 


Comp 

i 

T(l) 

T(2) 

Ordered  Failure  Times 

T(3)  T(4)  T(5) 

T(6) 

T(7) 

1 

186.409 

249.  124 

261.318 

362.  545 

406.506 

444. 204 

448. 009 

2 

59.693 

72.925 

75. 953 

88.763 

100. 177 

163.  702 

174.  788 

3 

43. 187 

81. 103 

83.922 

91.876 

93.052 

96.477 

98.  541 

4 

63. 788 

66. 730 

74. 303 

102.337 

102.  724 

104.502 

105.  736 

5 

25.015 

33. 172 

35.259 

35. 604 

57.  083 

61.971 

70.  905 

6 

12.920 

26. 710 

29.863 

29.998 

36.  680 

39.  804 

40. 992 

7 

15. 614 

23. 585 

26. 784 

31.  854 

33.  384 

35.  949 

37. 179 

8 

6. 164 

14. 714 

17.  345 

20. 198 

24. 394 

29.  221 

31.  666 

II.  Data  Summary 

Comp(i)  UT(i)  NC(i)  NF(i) 

TT(i) 

ELM(i) 

ER(i) 

ER(i)*Tr(l) 

1  10.0 

15 

7 

0.  18E+09 

0. 40E-07 

0. 23E-04 

4107. 8930 

2  10.0 

15 

7 

0.  43E+05 

0. 16E-03 

0.  96E-01 

4107.8930 

3  10.0 

15 

7 

0. 13E+11 

0.56E-09 

0. 33E-06 

4107. 8930 

4  10.0 

15 

7 

0. 69E+10 

0. lOE-08 

0. 60E-06 

4107.8930 

5  10.0 

15 

7 

0. 29E+05 

0. 24E-03 

0.  14E+00 

4107. 8930 

6  10.0 

15 

7 

0. 25E+06 

0. 28E-04 

0. 17E-01 

4107.8930 

7  10.0 

15 

7 

0. :5E+07 

0. 45E-05 

0. 27E-02 

4107. 8930 

8  10.0 

15 

7 

0. 41E+04 

0. 17E-02 

0.  lOE+01 

4107.8930 

III.  Estimation  Procedure  for  RSLOW 

Parameter 

df 

Value 

F.S 

0.92164 

ALPHA 

0. 1 

NFC 

56 

GHISQD 

112 

131.56 

(from  the  tabl*»  Chi 

-square  distribution  ) 

LMU 

0.00200 

RSLOW 

0.91712 

IV.  Workarea 

Comp  Ordered  Failure  Times  (h) 

raised  to 

the  power  of  Beta 

i  T’(l) 

T'(2)  T'(3) 

T‘(4)  T'(5)  T'(6)  T'(7) 

1 

0. 14E+07 

0. 30E+07 

0.  35E+07 

0.  84E+07 

0.  llE+08 

0.  14E+08 

0. 

15E+08 

2 

0. 68E+03 

0.  94E+03 

0.  lOE+04 

0. 13E+04 

0.  16E+04 

0.  34E+04 

0. 

38E+04 

3 

0.  25E+08 

0. 42E+09 

0. 49E+09 

0.  74E+09 

0.  79E+09 

0.  93E+09 

0. 

lOE+10 

4 

0.  64E+08 

0.  77E+08 

0.  12E+09 

0.  49E+09 

0.  50E+09 

0.  54E+09 

0.  57E+09 

5 

0.  37E+03 

0.  62E+03 

0.  69E+03 

0.  71E+03 

0.  17E+04 

0.  20E+04 

0. 

25E+04 

6 

0.  94E+03 

0. 66E+04 

0.  89E+04 

0.  90E+04 

0.  15E+05 

0.  19E+05 

0. 

21E+05 

7 

0.  77E+04 

0. 30E+05 

0.  45E+05 

0.  79E+05 

0.  92E+05 

0.  12E+06 

0. 

13E+06 

52 


8  0. 22E+02  0. 96E+02  0. 13E+03  0. 16E+03  0. 23E+03  0. 315+03  0. 35E+03 


Case  5  :  8  Weibull  components  in  parallel. 

- -  TEST  PLAN  1  -  Test  15  until  7  fails  for  each  component 

I.  Raw  Data 


Comp 

i 

T(l) 

TC2) 

Ordered  Failure  Times 

T(3)  T(4)  T(5) 

T(6) 

T(7) 

1 

3.  728 

4.  982 

5.226 

7.251 

8.  130 

8.  884 

8.960 

2 

2.296 

2.  805 

2.921 

3.414 

3.853 

6.  296 

6.  723 

3 

2.399 

4.506 

4.  662 

5. 104 

5.  170 

5.360 

5.475 

4 

4.  556 

4.  766 

5.307 

7.310 

7.337 

7.464 

7.553 

5 

2.  156 

2.860 

3.  040 

3.069 

4.  921 

5.  342 

6.  113 

6 

1.292 

2.671 

2.986 

3.  COO 

3.  668 

3.980 

4.  099 

7 

1.  763 

2.  663 

3.024 

3.596 

3.  769 

4.  059 

4.  198 

8 

0.  771 

1.  839 

2.  168 

2.525 

3.  049 

3.653 

3.958 

II.  Data  Summary 


CompCi)  UT(i) 

NCCi) 

NF(i) 

TT(i) 

ELM(i) 

ER(i) 

ER(i)*TT(i) 

1 

10.0 

15 

7 

0.  45E+04 

0.  16E-02 

0. 27E-01 

120.5522 

2 

10.0 

15 

7 

0.  24E+03 

0.  30E-01 

0.51E+00 

120.5522 

3 

10.0 

15 

7 

0.  27E+05 

0.  26E-03 

0. 45E-02 

120.5522 

4 

10.0 

15 

7 

0.  76E+05 

0.  92E-04 

0. 16E-02 

120.5522 

5 

10.0 

15 

7 

0.  32E+03 

0.  22E-01 

0. 38E+00 

120.5522 

6 

10.0 

15 

7 

0. 52E+03 

0.  13E-01 

0. 23E+00 

120.5522 

7 

10.0 

15 

7 

0. 13E+04 

0.  55E-02 

0.  95E-01 

120.5522 

8 

10.0 

15 

7 

0. 12E+03 

0.58E-01 

0. lOE+01 

120.5522 

III.  Estimation  Procedure  for  RSLOW 


Parameter  df  Value 


RS  0. 92659 

ALPHA  0.  1 

NFC  56 

CHISQD  116  135.89  (from  the  table  Chi-square  distribution  ) 

LMU  0. 07046 


RSLOW  0. 92598 


IV.  Workarea 


Comp  Ordered  Failure  Times  (h)  raised  to  the  power  of  Beta 
i  T'(l)  T'(2)  T'(3)  T'(4)  T'(5)  T'(6)  T'(7) 
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1 

2 

3 

4 

5 

6 

7 

8 


0.  35E+02 
0.  38E+01 
0.  52E+02 
0.  70E+03 
0. 41E+01 
0. 20E+01 
0. 63E+01 
0. 64E+00 


0. 77E+02 
0.52E+01 
0. 90E+03 
0. 86E+03 
0. 69E+01 
0. 14E+02 
0. 24E+02 
0. 28E+01 


0. 88E+02 
0.55E+01 
0. lOE+04 
0. 14E+04 
0.77E+01 
0. 19E+02 
0. 37E+02 
0. 37E+01 


0. 21E+03 
0.71E+01 
0. 16E+04 
0. 54E+04 
0. 78E+01 
0. 19E+02 
0. 65E+02 
0.48E+01 


0. 29E+03 
0. 86E+01 
0. 17E+04 
0. 55E+04 
0. 19E+02 
0. 32E+02 
0. 75E+02 
0. 66E+01 


0. 37E+03 
0. 19E+02 
0. 20E+04 
0. 60E+04 
0. 22E+02 
0. 40E+02 
0. 96E+02 
0. 90E+01 


0. 38E+03 
0. 21E+02 
0. 22E+04 
0. 63E+04 
0. 28E+02 
0. 44E+02 
0. llE+03 
0. lOE+02 


Case  8  :  8  Weibull  components  in  series-parallel. 

- -  (4  components  in  parallel  coimected  to  4  components  in  series) 

- -  TEST  PLAN  1  -  Test  15  until  7  fails  for  each  component 

I.  Raw  Data 

Comp 

Ordered  Failure  Times 

i 

T(l) 

T(2) 

T(3)  T(4) 

T(5) 

T(6) 

T(7) 

1 

5. 142 

6.  872 

7.209  10.001 

11.214 

12. 254 

12. 359 

2 

3. 184 

3.889 

4.  051  4.  734 

5.  343 

8.731 

9.322 

■  3 

3.344 

6.279 

6.497  7.113 

7.204 

7.469 

7.629 

4 

6.379 

6.673 

7.430  10.234 

10. 272 

10.450 

10.574 

5 

83. 383 

110.574 

117.532  118.681 

190.276 

206.571  236.351 

6 

25. 840 

53.419 

59.727  59.996 

73. 360 

79. 608 

81.983 

7 

24. 289 

36. 688 

41.663  49.551 

51.931 

55.920 

57.835 

8 

8.219 

19. 619 

23. 126  26. 931 

32.526 

38. 962 

42.221 

II.  Data  Summary 

Comp(i)  UT(i) 

NC(i) 

NF(i)  TTCi) 

ELM(i) 

ER(i) 

ER(i)*TT(i) 

1 

10.0 

15 

7  0. llE+05 

0.66E-03 

0. 38E-01 

399.5605 

2 

10.0 

15 

7  0. 40E+03 

0.  18E-01 

0. lOE+01 

399.5605 

3 

10.  0 

15 

7  0. 12E+06 

0.58E-04 

0. 33E-02 

399.5605 

4 

10.0 

15 

7  0. 33E+06 

0.21E-04 

0. 12E-02 

399.5605 

5 

10.  0 

15 

7  0. 26E+06 

0. 27E-C4 

0. 15E-02 

399.5605 

6 

10.0 

15 

7  0. 16E+07 

0. 44E-05 

0. 25E-03 

399.5605 

7 

10.0 

15 

7  0. 65E+07 

0. llE-05 

0. 61E-04 

399.5605 

8 

10.0 

15 

7  0. 67E+04 

0. lOE-02 

0. 60E-01 

399.5605 

III.  Estimation  Procedure  for  RSLOW 


Parameter  df  Value 


RS 

0.  93682 

ALPHA 

0. 1 

NFC 

56 

CHISQD 

104  122. 86 

LMU 

0.01922 

(from  the  table  Chi-square  distribution  ) 
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RSLOW 


0. 93382 


IV.  Workarea 


Comp 

i 

Ordered  Failure  Times  (h)  raised 
T'(l)  T'(2)  T'(3)  T'(4) 

to  the  power  of  Beta 
T'(5)  T'(6) 

T'(7) 

1 

0. 84E+02 

0. 18E+03 

0.  21E+03 

0.  51E+03 

0. 69E+03 

0. 88E+03 

0. 90E+03 

2 

0. 63E+01 

0. 87E+01 

0.  93E+01 

0. 12E+02 

0. 14E+02 

0. 32E+02 

0. 35E+02 

3 

0. 23E+03 

0. 40E+04 

0.  47E+04 

0. 71E+04 

0. 75E+04 

0. 88E+04 

0. 97E+04 

4 

0. 30E+04 

0.  37E+04 

0.  58E+04 

0.  23E-<-05 

0.  24E+05 

0.  26E+05 

0. 27E+05 

5 

0.  34E+04 

0.  57E+04 

0.  63E+04 

0.  65E+04 

0. 15E+05 

0. 18E+05 

0. 23E+05 

6 

0. 60E+04 

0.  42E+05 

0.  57E+05 

0.  58E+05 

0. 99E+05 

0.  12E+06 

0. 13E+06 

7 

0. 33E+05 

0.  13E+06 

0.  19E+06 

0.  33E+06 

0. 39E+06 

0.  49E+06 

0.55E+06 

8 

0. 36E+02 

0.  16E+03 

0.  21E+03 

0.  27E+03 

0. 37E+03 

0. 50E+03 

0.57E+03 
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APPENDIX  A.  DERIVATION  OF  FORMULA  FOR  BIAS  VALUE 

Suppose  Y  has  an  extreme-value  distribution;  i.e.,  Y ~  EV{a,  n). ' 
Suppose  ^  ^  ^(r+3)  ^  ^  the  middle  n-s  ordered 

statistics  in  a  random  sample  of  size  n  from  extreme-value  distribution 
with  pdf.  and  Cdf.  given  by 

J[y;  ii,a)  =  —  exp{  -  od  <  j;  <  c5o 

—  OQ  <  H  <  OO 

cr  >  0 

F(y;  iu,  ct)  =  1  -  exp{  (A. 2) 

where  cr  =  l/j? ,  —  In  6  and  In  is  the  natural  logarithm. 

Letting  denote  the  approximated  likelihood  function  and 
~  ix)/a,  from  Chapter  II  the  approximated  likelihood  equations 
for  fi  and  u  are  : 


AlnL  ^_5lnZ/ 
6/1  6/1 


(-4.3) 


5(1  -  %-s  +  Pn-s^in-s))  (°^i  “  ® 
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S  In  X/  ^  In  L  1  r  ^  I  V  /  Si  V  ^ 

=  "  "5"  “  ^^(r+1))  ~ 


6a 


6a 


{A  A) 


n—s 


5)(^  ^n—s  ^n— s'^(m— s))  Pi^{i)^^  ~  ^• 


It  is  not  possible  to  compute  the  conditional  s-bias  of  a  exactly  nor 
the  conditional  s-bias  of  ju.  The  conditional  s-bias  of  a  however,  can  be 
approximated  by  : 


rf  -iilLiill 

6a 


EL 


iiliLLli 

C  2  ^ 

6a 


{AS) 


Working  out  this  approximation  is  tedious,  but  Balakrishnan  and 
Varadan  [[Ref.  6  p-149]  provide  a  table  of  constants  for  the  conditional 
s-bias  of  a. 

Let  C  denote  the  bias  value  of  a  from  the  table.  Then 

EVal  -  g  ^ 
a  “  ^ 

E  Ca]  =  (T  -b  (jC  =  (1  -f  C)a. 

Thus  a*  =  al{\  -\r  C)  is  nearly  unbiased  estimator  of  a. 
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APPENDIX  B.  DERIVATION  OF  EQUATION  FOR  DEGREES  OF 

FREEDOM 

One  basic  expression  used  to  determine  the  degrees  of  freedom,  DF, 
for  the  chi-square  distribution  when  obtaining  confidence  limits  on  the 
reliability  of  a  complex  system  is  DF  =  2NFC,  where  NFC  is  total 
number  of  failed  components. 

In  chapter  V,  the  analysis  of  the  simulation  results  show  that  another 
formulae  for  computing  degrees  of  freedom  yielded  a  more  accurate  lower 
confidence  limit  procedure.  The  formula  DF  =  2NFC  -  NCOMP,  for 
example,  gfwe  more  accurate  results  than  procedures  that  used  DF  = 
2NFC.  Moreover,  the  accuracy  associated  with  any  fixed  formula  for  DF 
degraded  as  the  extent  of  truncation  increased.  That  is,  if  for  a  particular 
component,  NC  items  are  tesied  until  NF  fail,  the  accuracy  of  each  pro¬ 
cedure  decreased  as  NF  decreased.  The  decrease  in  accuracy  became  sig¬ 
nificant  for  values  .  ^  NF  <  NC/2.  Consequently,  a  formula  for  DF  was 
developed  that  included  NC/NF  as  one  of  the  terms  in  the  expression  for 
DF. 

This  formula  for  DF  has  the  form  DF  =  2NFC  +  c(NC/NF)  where 
NF  is  the  smallest  failure  truncation  value  for  each  component  type,  NC 
is  the  number  of  items  placed  on  test  for  that  component  with  smallest 
NF  and  c  is  an  unknown  constant.  An  equation  was  established  using  this 
formula  for  each  of  the  series  system  cases  simulated  using  the  value  of 
DF  that  yielded  the  most  accurate  result  and  the  appropriate  values  of 
NFC,  NC  and  NF  for  that  case.  Each  equation  can  be  solved  for  c.  Since 
there  were  thirty  series  system  cases  simulated  (See  Table  lA  through 
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Table  3B),  the  resulting  thirty  equations  yielded  thirty  values  of  c.  An 
averaging  process  using  least  square  methods  was  used  to  determine  one 
value  for  c.  The  following  set  of  equations  show  these  computations  for 
each  confidence  level  value. 

A.  ALPHA  =  0.1  (CONFIDENCE  90  %  ) 

152  =  160  +  c  reduces  to  -8  =  c 

232  =  240  +  c  reduces  to  -8  =  c 

168  =  176  +  1.36c  t-educes  to  -8  =  1.36c 

104  =  112  +  2.14c  reduces  to  -8  =  2.14c 

40  =  48  +  5.0c  reduces  to  -8  =  5.0c 
152  =  160  +  c  reduces  to  -8  =  c 

232  =  240  +  c  reduces  to  -8  =  c 

168  =  176  +  1.36c  reduces  to  -8  =  1.36c 

104  =  112  +  2.14c  reduces  to  -8  ==  2.14c 

40  =  48  +  5.0c  reduces  to  -8  =  5.0c 
144  =  160  +  c  reduces  to  -16  =  c 

224  =  240  +  c  reduces  to  -16  =  c 

160  =  176  +  1.36c  reduces  to  -16  =  1.36c 

96  =  112  +  2.14c  reduces  to  -16  =  2.14c 

64  =  48  +  5.0c  reduces  to  16  =  5.0c 
144  =  160  +  c  reduces  to  -16  =  c 

224  =  240  +  c  reduces  to  -16  =  c 

160  =  176  +  1.36c  reduces  to  -16  =  1.36c 

112  =  112  +  2.14c  reduces  to  0  =  2.14c 

64  =  48  +  5.0c  reduces  to  16  =  5.0c 
160  =  160  +  c  reduces  to  0  =  c 

256  =  240  +  c  reduces  to  16  =  c 


59 


192  =  176  +  1.36c  reduces  to  16  =  1.36c 

128  =  112  +  2.14c  reduces  to  16  =  2.14c 

64  =  48  +  5.0c  reduces  to  1 6  =  5.0c 
152  =  160  +  c  reduces  to  -8  =  c 

232  =  240  +  c  reduces  to  -8  ==  c 

176  =  176  +  1.36c  reduces  to  0  =  1.36c 
li2  =  112  +  2.14c  reduces  to  0  =  2.14c 
48  =  48  +  5.0c  reduces  to  0  =  5.0c 

The  above  equations  form  a  linear  system  that  is  overdetermined.  By 
using  the  normal  equation,  A^Ax  =  A'^b ,  we  then  find  the  least  square 
solution  to  that  overdetermined  system  which  yields  c  =  1.1.  Replacing 
the  constants  in  the  general  formula,  we  have  the  approximating  equation 
for  obtaining  the  degrees  of  freedom  ; 

DF  =  2NFC  +  1.1(NC/NF) 

B.  ALPHA  =  0.2  (CONFIDENCE  80  %  ) 

152  =  160  +  c  reduces  to  *8  =  c 

232  =  240  +  c  reduces  to  -8  =  c 

168  =  176  +  1.36c  reduces  to  -8  =  1.36c 

104  =  112  +  2.14c  reduces  to  -8  =  2.14c 

40  =  48  +  5.0c  reduces  to  -8  =  5.0c 
152  =  160  +  c  reduces  to  -8  =  c 

232  =  240  +  c  reduces  to  -8  =  c 

168  =  176  +  1.36c  reduces  to  -8  =  1.36c 

104  =  112  +  2.14c  reduces  to  -8  =  2.14c 

"  48  +  5.0c  reduces  to  -8  =  5.0c 
-  160  +  c  reduces  to  -16  =  c 

224  =  240  +  c  reduces  to  -16  =  c 
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160  =  176  +  1.36c  reduces  to  -16  =  1.36c 
96  -  112-1-  2.14c  reduces  to  -16  =  2.14c 
64  =  48  -f  5.0c  reduces  to  16  =  5.0c 
144  =  160  +  c  reduces  to  -16  =  c 

224  =  240  ■+  c  reduces  to  -16  =  c 

160  =  176  +  1.36c  reduces  to  -16  =  1.36c 
96  =  112  +  2.14c  reduces  to  -16  =  2.14c 
64  =  48  -1-  5.0c  reduces  to  16  =  5.Uc 
152  =  160  -F  c  reduces  to  -8  =  c 

240  ==  240  -1-  c  reduces  to  0  =  c 

176  =  176  +  1.36c  reduces  to  0  =  1.36c 
112  =  112-1-  2.14c  reduces  to  0  =  2.14c 
48  =  48-1-  5.0c  reduces  to  0  =  5.0c 
152  =  160  -1-  c  reduces  to  -8  =  c 

232  =  240  +  c  reduces  to  -8  =  c 

168  =  176  -1-  1.36c  reduces  to  -8  =  1.36c 
112  =  112  -1-  2.14c  reduces  to  0  =  2.14c 
48  =  48  -1-  5.0c  reduces  to  0  =  5.0c 

Using  the  normal  equation,  A''^A  x  =  A^b  again,  the  least  squares 
solution  is  c  =  1.0  .  The  approximate  equation  for  obtaining  the  degrees 
of  freedom  is  DF  =  2NFC  -1-  (NC/NF) 
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APPENDIX  C.  USERS'  GUIDE  FOR  RETP 


Reliability  Estimation  Test  Plan  (RETP). 


1.  Brief  Description 

RETP  is  a  computer  program  written  in  FORTRAN  that  runs 
on  the  Amdahl  mainframe  at  NPGS.  It  allows  the  user  to  simulate 
exponential  Weibull  failure  times  of  component  items  being  tested 
to  evaluate  the  accuracy  of  a  confidence  limit  estimation  procedure 
based  on  Type  II  data  censoring  (that  is,  testing  items  of  compo¬ 
nent  i  until  fi  of  them  fail). 

2.  Program  Input.  (INI.DATA) 

The  inputs  of  the  program  are  specified  to  the  program  via  an 
input  file  called  INI.DAT.  A  sample  input  file  is  shown  below. 

This  sample  input  refers  to  a  system  with  configuration  as  re¬ 
presented  by  Figure  4.4. 


This  file  contains  the  inputs  required  by  the  RETP  model. 

Update  only  the  numerical  values  between  dotted  lines  as  appropriate. 
Do  not  delete  any  of  the  comment  lines.  (INI.D/'"' 


Value 

Type 

Units 

Description 

Variable 

3 

INT 

Configuration  of  the  system 

1  =  all  subsystems  in  series 

2  =  all  subsystems  in  parallel 

3  =  subsystems  are  connected 

in  series-parallel 

CONFIG 

16807. 0 

REAL 

- 

Initial  random  seed 

ISEED 

15 

INT 

- 

Total  #  of  subsystems  in  system 

NCOMP 

10 

INT 

- 

#  of  exponential  subsystems 

NEXP 

5 

INT 

- 

#  of  Weibull  subsystems 

NWEI 

0 

INT 

- 

#  of  geometric  subsystems 

NGEO 

0. 10 

REAL 

- 

Desired  significance  level 

ALPHA 

1000 

INT 

- 

#  of  replications  desired 

NREP 

3 

INT 

Test  case  number 

1  =  all  exponential 

2  =  ell  Weibull 

3  =  EXP  +  WEI 

4  =  EXP  +  WEI  +  GEO 

TCN 

15 

INT 

• 

Number  of  cut  sets 

NCS 
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c  TEST  PLAN  :  Testing  NC(I)  items  of  component  i  until  NF(i)  of  them 

c  fail. 

c  (Use  REAL  numbers  ONLY  ! ! !  ) 

- 

c  Components  Configuration 

c  CSERIESCi)  CSSGROUP(i)  CSGKOUP(i)  NCGROUP(i) 

c - - - - — - - 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  8  14 

2  4  6  14 

2  4  6  14 

2  4  6  14 

2  4  6  14 

11  6  14 

11  6  14 

11  1  1 


c 


c  Corap  Comp  Comp  Parameters  Util  Test  Plan 

c  Number  Type  Scale  Shape  Time/Cycle  input 

c  I  TY(I)  PARM(1,1)  PARM(2,I)  UT(I)  UC(I)  NC(I)  NF(I) 


^  Int 

Int 

Real 

Real 

(hrs) 

Int 

Int 

Int 

1.0 

2.0 

0. 07500 

1.5 

10.0 

20.0 

15.0 

2.0 

2.0 

0. 05500 

2.0 

8.  0 

20.0 

15.0 

3.0 

2.0 

0. 17500 

1.5 

8.  0 

20.  0 

IS.  0 

4.0 

2.0 

O.C9500 

3.0 

2.0 

20.0 

15.0 

5.0 

2.0 

0. 12500 

2.0 

5.0 

20.0 

15.0 

6.  0 

1.0 

0. 05500 

1.0 

1,0 

20.0 

15.0 

7.0 

1.0 

0. 15000 

1.0 

1.0 

20.0 

15.0 

8.  0 

1.  0 

0. 15500 

1.0 

2.0 

20.  0 

15.0 

9.0 

1.0 

0.07500 

1.0 

2.  0 

20,0 

15.0 

10.0 

1.  0 

0. 05000 

1.0 

2.0 

20.  0 

15.0 

11.0 

1.0 

0. 17500 

1.0 

2.  0 

20.0 

15.0 

12.0 

1.0 

0. 08500 

1.0 

1.5 

20.0 

15.0 

13.0 

1.  0 

0. 17500 

1.0 

1.  0 

20.0 

15.0 

14.0 

1.  0 

0. 05000 

1.0 

A  *. 

20.0 

15.0 

15.0 

1.0 

0. 02500 

1.0 

1.0 

20.0 

15.0 

Note  ; 

TY'I)=1 

EXPONENTIAL 

P( surv) 

=  exp(-PARM(l,I) 

*  T) 

c  TY(I)=2  WEIBULL  P(surv)  =  exp('(PARM(l,I)*T)**PARM(2,i)) 

c  Ty(I)=3  GEOMETRIC  P(surv)  =  PARM(1,I)  **  T 

c - - - - - - 

c  SYSTEM  CONFIGURATION  :  Identification  of  CUT  SETS 

c  Minimum  groups  of  subsystems  that  have  to  fail 

c  for  the  system  to  fall. 

- - - - - - - - - 

c  Cut  Set  #  in  Set  List  of  Components  in  Cutset 
c  J  C0MP(J,1)  C0MP(J,2) . up  to  C0MP(J,1)  components 

c - - - - - 
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1 

1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


10  0 
2  0  0 

3  0  0 

1  4  0  0 

1  5  0  0 

1  6  0  0 

1  7  0  0 

1  8  0  0 

1  9  0  0 

1  10  0  0 

1  11  0  0 

1  12  0  0 

1  13  0  0 

1  14  0  0 

1  15  0  0 


■)  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


3.  Program  Flow  and  Logic.  (NAME.DEF,PARM.DEF  and 
RETP.FOR) 

Input  parameters  are  first  read  in  by  the  program  by  calling  the 
INPUT  subroutine.  The  program  then  evokes  the  SIM  subroutine 
which  generates  the  random  failure  times  and  computes  the  key  sta¬ 
tistics  required  in  the  procedure.  The  next  subroutine  EVAL  deter¬ 
mines  the  measures  of  accuracy  for  each  case.  REPORT  is  the 
subroutine  which  generates  the  output  file  for  the  run  OUT.DAT. 

The  variables  in  the  program  RETP.FOR  are  described  in  the  file 
NAME.DEF  as  listed  below. 


0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


This  file  contains  the  declaration  for  input  and  output  variables 
used  in  the  RETP  model.  (NAME.DEF) 


Input  Variables. 


ISEED  =  initial  random  seed  selected.  (Real) 

SEED  =  current  random  seed  (Real) 

RS  =  true  overall  system  reliability. 

ALPHA  =  desired  significance  level  (Real) 

NREP  =  number  of  replication  desired  for  simulation  (Integer) 

TPN  =  test  plan  number 

TCN  =  test  case  number  (1,  2,  3  or  4) 

NCOMP  =  total  #  of  components  in  the  system  (Integer) 

I^EXP  =  number  of  components  with  EXP  failure  times.  (Integer) 

NWEI  =  number  of  components  with  WEI  failure  times.  (Integer) 

NGEO  =  number  of  components  with  GEO  failure  times.  (Integer) 

Distribution:  Exponential  WEibull  GEOmetric 

Ty(i)  =  type:  1  2  3 

PARM(l,i)  :  Scale(l/hr)  Scale(l/hr)  Prob 

PARM(2,i)  :  -  Shape 

UT(i)  =  utilization  time  (hrs)  for  component  i  (EXP  and  WEI) 

UC(i)  =  utilization  cycles  for  component  i  (GEO  only) 

NC(i)  =  number  of  test  samples  (sample  size)  for  component  i. 

NF(i)  =  desired  niunber  of  failures  in  test  for  component  i. 
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c  NCS  =  number  of  cut  sets  for  the  system. 

c  C0MP(J,K)  =  kth  parameter  of  cutset  j  (  first  being  the  number  of 
c  components  belonging  to  the  cut-set) 

c  CONFIG  =  configuration  number  of  components  arrangement( integer), 
c  CSERIES(i)  =  number  of  subsystems  in  series  in  a  subsubgroupC integer) 
c  CSGROUP(i)  =  number  of  subsystems  i  in  a  subgroup  (integer) 
c  CSSGROUP(i)=  number  of  subsystems  i  in  a  subsubgroup  (integer) 
c  NCGROUP(i)  =  number  of  subsystems  i  in  a  group  (integer) 
c 

c  Assumed  Variables. 

c  - - 

c  MAXCOMP  =  maximum  number  of  components  allowed  in  the  system 
c  MAXREP  =  maximum  number  of  replication  permitted, 
c  MAXCUT  =  maximum  number  of  cut-sets, 
c 

c  Program  and  Output  Variables. 

c  - - - - - 

c  TT(i)  =  total  accumulated  failure  time  (in  hour)  for  component  i 
c  (EXP  and  WEI  only) 

c  TC(i)  =  total  accumulated  cycles  to  failure  (incl.  failure  cycle) 
c  for  component  i  (GEO  only) 

c  EEETA(i)  =  estimate  for  shape  parameter  of  component  i  (if  Weibull) 
c  RELl(j)  =  actual  reliability  for  cut-set  j. 

c  REL2(j)  =  computed  reliability  for  cut-set  j  for  current  replication 

c  ELM(i)  =  estimated  component  failure  rate  (I'hrs)  for  component  i 

c  ELMAX(m)  =  max.  estimated  component  failure  rate  for  rep.  m  (1/hrs) 

c  ER(i)  =  ratio  of  estimated  failure  rate  to  ELMAX. 

c  NFC(m)  =  total  number  of  failed  test  components, 

c  LhU(m)  =  upper  confidence  limit  for  failure  rate  (1/hrs). 

c  RSL(m)  =  lower  confidence  limit  estimated  for  system  reliability 

c  for  the  mth  replication, 

c  ORSL(m)  =  ascending  order  of  RSL(M). 
c  RSLOW  =  (l-ALPHA)xlOO  percentile  of  set  of  RSL(M). 
c  LEVEL  =  achieved  confidence  interval  i. e  proportion  of  RSL(M) 
c  that  are  lesser  than  RS  (cons'trvative  estimate) 

c 

c—  END  OF  NAME.DEF . 

c 

Together  with  the  main  program  in  RETP.FOR  are  the  other  sub¬ 
routines  needed  in  the  simulation.  The  declaration  of  variables  is  done  in 
the  file  PARM.DEF.  Relevant  descriptions  are  included  as  comment  lines 
in  the  source  code  to  '-^Ip  explain  the  program  segments.  A  listing  of 


PARM.DEF  and  RETP.FOR  is  given  below. 


c  This  file  contains  the  declaration  for  input  and  output  variables 
c  used  in  the  RETP  model.  (PARM.DEF) 


c 

INTEGER  MAXCOMP,  MAXREP 

PARAMETER(  MAXCOMP  =  100,  MAXREP  =  1000,  MAXCUT  =  20  ) 
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REAL*8  I SEED,  SEED 

INTEGER  NREP,  TCN,  NCOMP,  NEXP,  NWEI,  NGEO,  NCS, 

+  NC(MAXCOMP),NF(MAXCOMP),  TY(MAXCOMP),  NFC(MAXREP), 

+  UC(KAXCOMP),  TC(MAXCOMP),  COMP(MAXCUT,MAXCOMP) , 

+  CONFIG,  CSERIES(MAXCOMP),  CSSGROUP(MAXCOMP) , 

+  CSGROUP(MAXCOMP) ,NCGROUP(MAXCOMP) ,DFR 

REAL*8  RS,  ALPHA,  UT(MAXCOMP),  TT(MAXCOMP), 

+  PARM(2,MAXCOMP),  ELM(MAXCOMP) ,  ER(MAXCOMP), 

+  LMU(MAXREP),  RSL(MAXREP),  ORSL(MAXREP) , 

+  ELMAX(MAXREP),  RSLOW,  LEVEL,  EBETA(MAXCOMP) , 

+  RELl(MAXCUT),  REL2(MAXCUT) 

c 

COMMON/BLOCKJ/ISEED,  SEED,  NREP,  TCN,  NCOMP,  NC,  NF,  NEXP,  NWEI, 
+  NGEO,  NCS,  TY,  NFC,  UC,  TC,  COMP,  CONFIG,  CSERIES, 

+  CSSGROUP,  CSGROUP,  NCGROUP,  DFR 

C0MM0N/BL0CK2/RS,  ALPHA,  UT,  TT,  PARM,  ELM,  ER,  LMU, 

+  RSL,  ORSL,  ELMAX,  RSLOW,  LEVEL,  EBETA,  RELl,  REL2 

c 

c  END  OF  PARM.DEF . - . . 

c 

c  This  file  contains  the  main  program  and  the  subroutines 
c-  for  the  Reliability  Estimation  Test  Plan  (RETP)  model, 
c 

c  IBM  Mainframe  version. 

c  Test  Plan  1  ;  Testing  NC(I)  items  for  component  i 
c 

C-- ---------- ---------- — 

c  1.  Main  Program  (RETP). 

PROGRAii  RETP 
c 

c  Include  the  declaration  files, 
c 

INCLUDE  'NAME  DE7’ 

INCLUDE  'PARM  DEF' 
c 

c  Read  the  input  data, 
c 

CALL  INPUT 
c 

c  Activate  simulation 
c 

CALL  SIM 
c 

c  Process  and  evaluate  output  data, 
c 

CALL  EVAL 
c 

c  Generate  simulation  report. 

G 

CALL  REPORT 
c 

STOP 

r'T) 
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c  2.  Input  Initialisation  Subroutine  (INPUT). 


SUBROUTINE  INPUT 
c 

c  This  subroutine  reads  in  the  inputs  for  the  RETP  model, 
c 

c  Include  the  declaration  file, 
c 

INCLUDE  'PARK  DEF' 
c 

INTEGER  I,  J,  K,  DUM2(11),  DUM3(4) 

REAL*8  DUM1(7) 
c 

c  Read  data  from  file  'INI. DAT*  designated  as  logic  unit  1. 
c 

OPEN  CUNIT=1,FILE=’/INI  data  Al') 
c 

READCl,*)  CONFIG 
READd,*)  ISEED 
READ(1,*)  NCOMP 
READd,*)  NEXP 
READCl,*)  NWEI 
READd,*)  NGEO 
READCl,*)  ALPHA 
READCl,*)  NREP 
READCl,*)  TCN 
READCl,*)  NCS 
c 

DO  40  I  =  1,  NCOMP 
READCl,*)  DUM3 
CSERIES(I)  =  DUM3(1) 

CSSGROUP(I)  =  DUM3(2) 

CSGROUP(I)  =  DUM3(3) 

NCGROUP(I)  =  DUM3(4) 

40  CONTINUE 
c 

DO  50  K  =  1,  NCOMP 
READCl,*)  DUMl 
c 

I  =  NINTCDUMld)) 

TYCI)  =  NINT(DUM1(2)) 
c 

IF  (TYCI).EQ.  1)  THEN 
PARM(1,I)  =  DUM1(3) 

PARM(2,I)  =  DUM1(4) 

UT(I)  =  DUM1(5) 

NCCI)  =  NINT(DUM1(6)) 

NF(I)  =  NINT(DUM1(7)) 

EBETA(I)  =  0 

ELSEIF  (TY(I).EQ.2)  THEN 
PARM(1,I)  =  DUM1(3) 

PARM(2,I)  =  DUM1(4) 

UTCI)  =  DUM1(5) 

NC(I)  =  NINT(DUM1C6)) 

NF(I)  =  NINT(DUM1(7)) 

ELSEIF  (TYCD.EQ.B)  THEN 
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PARM(1,I)  =  DUM1(3) 
PARM(2,I)  =  DUM1(4) 
UC(I)  =  NINT(DUiIl(5)) 
UT(I)  =  DUM1(5) 

NC(I)  =  NINT(DUM1(6)) 
NF(I)  =  NINT(DUM1(7)) 
ENDIF 
50  CONTINUE 

DO  80  I  =  1,  NCS 
READCl,*)  DUM2 
J  =  DUM2(1) 

C0MP(J,1)  =  DUM2(2) 

DO  70  K=1,C0MP(J,1) 

C0MP(J,K+1)  =  DUM2(K+2) 
70  CONTINUE 
80  CONTINUE 

CL0SECUNIT=1) 

RETURN 

END 


c  3.  Subroutine  for  Simulation. 


SUBROUTINE  SIM 
c 

c  This  subroutine  simulates  NPRP  possible  outcomes  of  the  test  plan 
c  desired  in  order  to  obtain  the  raw  estimates  of  IiMU(M)  and  RSL(M) 
c  for  each  of  the  replication, 
c 

c  Include  the  declaration  file 
c  and  declare  local  variables, 
c 

INCLUDE  'PARM  DEF' 
c 

INTEGER  I,  J,  K,  M,  ISUM,  KEY,  L,  C,  CC,  CCC 
REAL*8  UNI 

REAL*8  SUM,  PROD,  RP,  RPP,  RSP,  RSS, 

+  FT(MAXCOMP),  0FT(MAXC0MP),YY(MAXC0MP) 

c 

LOGICAL  TEMP 
c 

SEED  =  ISEED 
c 

c  Compute  overall  true  system  reliability  RS. 
c 

C  =  0 
CC  =  0 
CCC  =  0 
RP  =  0.  0 
RPP  =  1.0 
RSP  =  0. 0 
RSS=  1.0 
RS  =  1.  0 
J  =  1 

TEMP  -  .FALSE. 

DO  WHILE(.  NOT. TEMP) 
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DO  20  L  =  J,  J  +  (GSERIES(J)  -  1) 

PROD  =  1.0 

DO  10  I  =  1,  C0MP(L,1) 

K  =  C0MP(L,I+1) 

PROD  =  PR0D*(1-  SURV(TY(K),PARM(1,K),PARM(2,K),UT(K))) 
10  CONTINUE 

RELIC L)  =  1.  0  -  PROD 
RSS  =  RSS  *  RELl(L) 

20  CONTINUE 
L  =  L  -  1 

RP  =  1.0  -  (1.0  -RSS)*(1.  0  -RP) 

RSS  =1.  0 

IF  (L.GE.  (C  +  CSSGROUPCL)))  THEN 
RPP  =  RPP  *  RP 
RP  =  0.  0 

C  =  C  +  CSSGROUPCL) 

ENDIF 

IF  (L.GE. (CC  +  CSGROUPCL)))  THEN 

RSP  =  1.  0  -  (1.  0  -  RPP)*(1.  0  -  RSP) 

RPP  =  1.  0 

CC  =  CC  +  CSGROUP(L) 

ENDIF 

IF  (L.GE. (CCC  +  NCGROUP(L)))  THEN 
RS  =  RS  *  RSP 
RSP  =  0.  0 

CCC  =  CCC  +  NCGROUP(L) 

ENDIF 

J  =  J  +  CSERIES(L) 

IF  (L.GE.NCOMP)  TEMP  =  .TRUE. 

ENDDO 

c 

c  start  of  Simulation 
c  (Initialize  replication  counter  M) 
c 

M  =  1 

DO  WHILE  (M.LE.NREP) 

o 

c  Test  Plan  :  Sample  and  determine  imknown  TT(I) 
c  ------ —  with  known  NC(I)  until  NF(I)  fails. 

c 

c  Generate  NC(I)  failure  times,  put  them  in  ascending  order 
c  with  the  smallest  failure  time  on  the  top  of  the  list, 
c 

DO  70  I  =  1,  NCOMP 
c 

DO  40  K  =  1,  NC(I) 

CALL  LPJ®PC(SEED,UNI,1) 

IF  (TY(I).EQ. 1)  THEN 
FT(K)=  "L0G(UNI)/PARM(1,I) 

ELSEIF  (Ty(I).EQ.  2)  THEN 

FT(K)=( 1. 0/PARM( !,!))*( -LOG(UNI) )**( 1. 0/PARM( 2,1)) 
ELSEIF  CTY(I).EQ.3)  THEN 
FT(K)=  1.0 

DO  WHILE  (UNI.LT.PARMd,!)) 

FT(K)  =  FT(K)  +  1.  n 
CALL  LRNDPC(SEED,U:.x,l) 
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o  o 


ENDDO 

ENDIF 

40  CONTINUE 

c 

c  Buble  Sort  the  failure  times  in  ascending  order, 
c 

CALL  BUBBLE (NC( I), FT, OFT) 
c 

c  Take  logarithm  of  the  ordered  failure  times, 
c 

DO  45  K=l,  NC(I) 

YYCK)  =  LOGCOFTCK)) 

45  CONTINUE 
c 

c  Compute  the  total  time  accumulated  in  the  test  and  the  estimate 
G  for  the  failure  rate  of  the  component  as  in  the  procedure. 


IF  (TY(I).NE.2)  THEN 
SUM  =  0.  0 
DO  50  K  =  1,  NF(I) 

SUM  =  SUM  +  OFT(K) 

50  CONTINUE 

TT(I)  =  FL0AT(NC(I)'NF(I))*0FT(NF(I))  +  SUM 
IF  (TY(I).EQ.  1)  THEN 

ELM(I)  =  FL0AT(NF(I)  -  l)/Tr(I) 

ELSEIF  (TYCI).EQ.3)  THEN 

ELM(I)  =  FL0AT(NT(I)  -  l)/Tr(I) 

ENDIF 

ELSEIF  CTY(I).EQ. 2)  THEN 
c 

CALL  APR0XMLE(YY,NC(I) ,NF(I) ,EBETA(I)) 
c 

EBETA(I)  =  BN(NC(I))*EBETA(I) 
c 

SUM  =  0. 0 

BO  60  K  =  1,  NF(I) 

SUM  =  SUM  +  0FT(K)**i.BETA(I) 

60  CONTINUE 
c 

TT(I)  =  FLOAT(NC(I)-NF(I'>)*OFT(NF(I))**EBETA(I)  +  SUM 
ELM(I)  =  FL0AT(NF(I))/lTa) 

ENDIF 

c 

70  CONTINUE 
c 

c  Determine  the  total  number  of  failed  test  items, 
c 

ISUM  =  0 

DO  80  I  =  1,  NCOMP 
ISUM  =  ISUM  +  NF(I) 

80  CONTINUE 

NFC(M)  =  ISUM 
c 

c  Determine  the  maximiim  failure  rate  estimate 
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c  and  identify  the  component, 
c 

ELMAX(M)  =  0.0 
KEY  =  0 

DO  90  I  =  1,  NCOMP 

IF  (ELM(I).GT.ELMAX(M))  THEN 
ELMAX(M)  =  ELMCI) 

KEY  =  I 
ENDTF 
90  CONTINUE 
c 

c  Compute  the  ratios  of  the  failure  rate  estimate  to  their  maximum 
c 

DO  100  I  =  1,  NCOMP 

ER(I)  =  ELM(I)/ELMAX(M) 

100  CONTINUE 
c 

c  Determine  LMU(M) 
c 

SUM  =  0.  0 

DO  110  I  =  1,  NCOMP 

SUM  =  SUM  +  (ER(I)*Tr(I)) 

110  CONTINUE 

DFR  =  2*NFC(M)-2*NC0MP 
LMU(M)  =  CHISQDC 1-ALPHA, DFR)/C2*SUM) 

Compute  estimate  of  overall  reliability  RSL(M)  for  the  system. 

C  =  0 
CC  =  0 
CCC  =  0 
RP  =  0.  0 
RPP  =  1.  0 
RSP  =  0. 0 
RSS=  1.0 
RSLCM)  =  1.0 
J  =  1 

TEMP  =  .FALSE. 

DO  WHILE(. NOT. TEMP) 

DO  120  L  =  J,  J  +  (CSERIES(J)  -  1) 

PROD  =  1.  0 

DO  115  1=1,  COMP(L,l) 

K  =  COMP(L,I+l) 

IF  (TY(K).EQ.  1)  THEN 

PROD  =  PROD*(l-SURV(TY(K),LMU(M)*ER(K),EBETACK),UT(K))) 
ELSEIF  (TY(K).EQ.  2)  THEN 

PROD  =  PROD*(l-SURV(TY(K),(LMU(M)*ER(K))**(l./EBETA(K)), 

+  EBETA(K),UT(K))) 

ELSEIF  (TY(K).EQ.3)  THEN 

PROD=PROD*(l-SURV(TY(K) , 1.  DO-LMU(M)*ER(K) ,0.  D0,UT(K))) 
ENDIF 
115  CONTINUE 

REL2(L)  =  1.  0  -  PROD 
RSS  =  RSS*REL2(L) 

120  CONTINUE 
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L  =  L  -  1 

RP  =  1.  0  -  (1.0  -  RSS)*(1.  0  -  RP) 

RSS  =1.  0 

IF  (L.  GE.  (C  +  CSSGROUPCL)))  mN 
RPP  =  RPP  *  RP 
RP  =  0.0 

C  =  C  +  CSSGROUPCL) 

ENDIF 

IF  (L.GE. (CC  +  CSGROUP(L)))  THEN 

RSP  =  1.0  -  (1.0  -  RPP)  *  (1.0  -  RSP  ) 

RPP  =  1.  0 

CC  -  CC  +  CSGROUP(L) 

ENDIF 

IF  (L.GE.CCCC  +  NCGROUP(L)))  THEN 
RSL(M)  =  RSL(M)  *  RSP 
RSP  =  0.  0 

CCC  =  CCC  +  NCGROUP(L) 

ENDIF 

J  =  J  +  CSERIES(L) 

IF  (L.GE.NCOHP)  TEMP  =  .TRUE. 

ENDDO 

Increment  replication  counter, 
c 

M  =  M  +  1 
c 

ENDDO 

c 

RETURN 

END 

c----------- 

c  4.  Subroutine  for  Evaluation. 

c--- — - - - - — - -  —  ------  —  ---- 

SUBROUTINE  EVAL 
c 

c  This  subroutine  calls  BUBBLE  to  sort  the  array  RSL(NREP)  in 
c  ascending  order  to  get  an  ordered  array  ORSL(NREP).  It  also 
c  determine  the  estimate  for  RSLOW  at  the  specified  significance 
c  level  ALPHA  and  the  value  of  Level  in  which  ORSL(LEVEL)  is  closest 
c  to  the  true  reliability  RS. 
c 

c  Include  the  declaration  files 
c  and  declare  the  local  variables, 
c 

INCLUDE  'PARM  DEF' 
c 

INTEGER  INDEX,  M 
REAL*8  DIFF 
c 

c  Order  the  array  RSL(NREP)  in  ascending  order, 
c 

DO  10  M  =  1,  NREP 
ORSL(M)  =  RSL(M) 

10  CONTINUE 
c 
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c  Bubble  Sort.  Sink  the  larger  of  the  pair, 
c 

CALL  BUBBLE(NREP,RSL,ORSL) 
c 

c  Determine  the  (1-ALPHA)  %  lower  confidence  bound  for  the  system 
c  reliability, 
c 

RSLOW  =  0RSL(NINT(NREP*(1 -ALPHA))) 
c 

c  Finding  the  %  confidence  level  for  the  true  reliability  RS. 
c  (i. e  the  proportion  of  RSL(M)  lesser  than  RS) 
c 

DIFF  =  1.  0 
INDEX  =  0 

DO  200  M  =  1,  NREP 

IF  (ABS(0RSL(M)  -  RS).LT.DIFF)  THEN 
DIFF  =  ABS(0RSL(M)  -  RS) 

INDEX  =  M 
ENDIF 
200  CONTINUE 
c 

LEVEL  =  FLOATC INDEX) /NREP 
c 

c  Recoi.1  evaluated  parameters  in  RAW1.DAT  (unit  2). 
c 

0PEN(UNIT=2) 

WRITE(2  300) 

300  F0RMAT(ix,  '  M  LMU(M)  ELMAX(M)  RSL(M)*, 

+  '  ORSL(M)  NFC(M)’) 

DO  500  M  =  1,  NREP 

WRITE(2,400)  M,LMU(M),Ej  .xAX(M) ,RSL(M) ,0RSL(H) ,NFC(M) 

400  F0RMAT( lx, 16 , 2F12.  7 , 2F12.  7,16) 

500  CONTINUE 

CL0SE(UNIT-2) 

c 

RETURN 

END 

c---------- - - - - - - 

c  5.  Subroutine  for  Report  Generation. 

c - - - 

SUBROUTINE  REPORT 
c 

c  This  subroutine  record  the  simulation  result  into  the  'OUT  DATA' 
c  file  indicated  by  logic  3. 
c 

c  Include  the  declaration  files 
c  and  declare  local  variables, 
c 

INCLUDE  '  PARM  DEF' 

INTEGER  I,  J,  K,  DUM(IO) 
c 

c  Write  to  output  file  'OUT  DATA'  designated  as  logic  unit  3. 
c 

0PEN(UNIT=3) 

c 

WRITE(3,10) 
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WRITE(3,20)  NREP 
WRITE(3,25)  NCOMP, CONFIG 
IF  (CONFIG. EQ. 1)  THEN 
WRITE(3,26) 

ELSEIF  (CONFIG. EQ. 2)  THEN 
WRITE(3,27) 

ELSE 

WRITE(3,28) 

ENDIF 

WRITE(3,29)  DFR 
WRITE(3,30) 

WRITE(3,40) 

WRITE( 3,50)  I SEED , NCOMP .ALPHA, NCS ,TCN 

WRITE(3,60) 

DO  200  I  =  1,  NCOMP 

WRITE(3,70)  I,TY(I),PARM(1,I),PARM(2,I),UT(I),NC(I),NF(I) 

200  CONTINUE 
WRITE(3,80) 

WRITE(3,90) 

DO  300  I  =  1,  NCOMP 

WRITE(3 , 100)  I ,NF( I) ,TT( I) ,ELM( I) ,ER( I) ,EBETA( J) 

300  CONTINUE 

WRITE(3,110) 

WRITE(3,120) 

DO  500  J  =  1,  NCS 
DO  400  K  =  1,  10 
DUM(K)  s  COMP(J,K) 

400  CONTINUE 

WRITE(3,130)  J,DUM,REL1(J),EEL2(J) 

500  CONTINUE 

WRITE(3,140) 

WRITE(3,150)  RS,ELMAX(NREP),LMU(NREP),RSLOW, LEVEL 

10  FORMAT(lx,’OUTl  DATA  :  Output  File  of  the  RETPl  simulation') 

20  FORMATdx,'  after  ',15,'  replication',/) 

25  FORMATdx, 'COMMENTS  :  '  ,12, '  COMPONENTS  IN  CONFIGURATION:  '  ,12  ) 

26  FORMATdx,'  (SERIES  SYSTEM)  ',/) 

27  FORMATdx,'  (PARALLEL  SYSTEM)  ',/) 

28  FORMATdx,’  (SERIES-PARALLEL  SYSTEM)  ',/) 

25  FORMATdx,'  DF  =  2  *  (NFC  -  NCOMP)  =',I4  ) 

30  FORMAT( lx, 'Input  parameter: ' ,/) 

40  FORMATdx,'  I  SEED  NCOMP  ALPHA  NCS  TCN',/) 

50  F0RMAT(lx,F10. 1,18, F8. 4,216,/) 

60  FORMATdx,'  I  TY(I)  PARMl(I)  PARM2(I)  UT(I)  NC(I)  NF(I)’,/) 
70  FORMAT(lx,I2.I6,lx,2F9.5,F8.2,2x,2I6) 

80  FORMATdx,/, 'Output  Parameters  for  the  Last  Replication:',/) 

90  FORMATdx,’  I  NF(I)  TT(I)  ELM(I)  ER(I)  ', 

+  *  EBETA(I)',/) 

100  F0RMAT(lx,I2.I6,4x,E14.  7,2x,F9. 7,2x,F9.  7,2x,F9.  7) 

110  FORMATdx,/,  Tut-Set  Data:  ’ ,/) 

120  FORMATdx,'  J  NUM  Component  List 

+  '  RELl(J)  REL2(M)  ’ ,/) 

130  FORMAT(lx,I2,I5,9I3,2F12.9) 

140  FORMAT(lx,/,'  RS  ELMAX(M)  LMU(M)  ’ 

+  '  RSLOW  LEVEL  ',/) 
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o  o 


150  F0RMATClx,5F12. 7,/) 

CLOSECUNIT  =  3) 

RETURN 

END 


c  This  portion  of  the  file  contains  functions  and  subroutines 
c  used  in  the  RETP  model. 


c  A.  Random  Number  Generating  Subroutine  (LRNDPC). 


SUBROUTINE  LRNDPC  (DSEED,U,N) 

INTEGER  N,  I 
REAL*8  U(N) 

REAL*8  D31M1,  DSEED,  D31 


DATA  D31M1/2147483347.D0/ 

DATA  D31  /2147483648.D0/ 

DO  5  I  =  1,  N 

DSEED  =  DMODC 16807.  D0*DSEED,D31M1) 
5  U(I)  =  DSEED  /  D31 
RETURN 
END 


B.  Survivability  Function. 


FUNCTION  SURVCTYPE, PARI, PAR2, UTIL) 
c 

c  This  function  returns  the  survival  probability  of  the  component  of 
c  different  types  (TYPE)  with  scale  (PARI)  and  shape  (PAR2)  parameters 
c  given  the  specified  utilization  times  or  cycles  (UTIL), 
c 

INTEGER  TYPE,  N 
REAL*8  PARI,  PAR2,  UTIL 
c 

IF  (TYPE.EQ.  1)  THEN 

SURV  =  EXP(-(PAK1*UTIL)) 

ELSEIF  (TYPE.EQ.  2)  THEN 

SURV  =  EXP(-((PAR1*UTIL)**PAR2)) 

ELSE 

N  =  NINT(UTIL) 

SURV  =  PAR1**N 
ENDIF 
END 

- - - - - - — - - - - - - 

c  C.  Bubble  Sort  Routine  in  Ascending  Order. 

Q - - - - - - - - - 

SUBROUTINE  BUBBLE (N, LI ST, OLI ST) 
c 

c  This  subroutine  performs  a  bubble  sort  in  increasing  order, 
c  (i.  e  sink  the  greater  numeral)  for  the  first  N  terms  in  an  array 
c  LIST  and  returns  the  result  in  OLIST. 
c 
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LOGICAL  DONE 
INTEGER  N,  K,  PAIR 
REAL  TEMP 

REAL*8  LISTC*),  OLIST(*) 


c  Sink  the  larger  of  the  pair, 
c 

DO  50  K  -  1,  N 

OLIST(K)  =  LIST(K) 

50  CONTINUE 

PAIR  =  N  "  1 
DONE  =  .FALSE. 

DO  WHILE  (.NOT. DONE) 

DONE  =  .TRUE. 

DO  100  K  =  1,  PAIR 

TF  (OLIST(K).GT.OLIST(K+l))  THEN 
TEMP  =  OLIST(K) 

OLIST(K)  =  OLISTCK+ll 
0LIST(K+1)  =  TEMP 
DONE  =  .FALSE. 

ENDIF 

100  CONTINUE 

PAIR  =  PaIR  -  1 
ENDDO 
RETURN 
END 


c  D.  Unbiasing  Factor  for  Biased  MLE  for  Weibull  Shape  Parameter. 


FUNCTION  BN(I) 
c 

c  This  function  returns  the  value  of  the  unbiased  factor  for  the  biased 
c  approximate  MLE  of  the  shape  parameter  of  a  Weibull  distribution 
c  with  a  sample  size  of  N. 


INTEGEi’  I 

IF  (I.LE.5)  TliEN 

BN  =  (1*0. 699)/(5. 0) 

ELSEIF  (I.EQ.6)  THEN 
BN  =  0.  752 

ELSEIF  (I.EQ.  7)  THEN 
BN  =  0.  786 

ELSi  .:I.EQ.  8),  THEN 
BN  =  0.  82 

ELSEIF  (I.EQ.  9)  THEN 
BN  =  0. 8395 

FLSEIF  (I.EQ.  10)  THEN 
BN  =  0.  859 

ELSEIF  (I.EQ.  11)  THEN 
BN  =  0.871 

ELSEIF  (I.EQ. 12)  THEN 
BN  =  0.  883 
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ELSEIF  (I.EQ.  13)  THEN 
BN  =  0.  892 

ELSEIF  (I.EQ. 14)  THEN 
BN  =  0.  901 

ELSEIF  (I.EQ.  15)  THEN 
BN  =  0. 9075 
ELSEIF  (I.EQ. 16)  THEN 
BN  =  0.914 

ELSEIF  (I.EQ.  17)  THEN 
BN  =  0.  9185 
ELSEIF  (I.EQ.  18)  THEN 
BN  =  0.  923 

ELSEIF  (I.EQ.  19)  THEN 
BN  =  0.  927 

ELSEIF  (I.EQ.  20)  THEN 
BN  =  0.  931 

ELSEIF  (I.LE. 25)  THEN 

BN  =  0.931  +  (I-20)-V0.  014/5.  0 
ELSEIF  (I.LE.  30)  THEN 

BN  =  0.945  +  (I-25)*0.  01/5.0 
ELSEIF  (I.LE.  40)  THEN 

3N  =  0.955  +  (I-30)*0. 011/10. 0 
ELSEIF  (I.LE.  60)  THEN 

BN  -  0.966  +  (I-40)*0.  012/20.  0 
ELSEIF  (I.LE. 80)  THEN 

BN  =  0.978  +  (I-60)*0. 006/20. C 
ELSEIF  (I.LE,  100)  THEN 

BN  =  0.984  +  (I-80)*0.  003/20. 0 
ELSEIF  (I.LE.  120)  THEN 

BN  =  0.987  +  (I-*100)*0.  003/20.0 
ELSE 

BN  =  1.  0 
ENDIF 
RETURN 
END 


c  E. Unbiasing  Factor  for  approx  MLE  for  scale  parameter. 


FUNCTION  BIAS(N,S) 
c 

c  This  fiinction  returns  th^-  value  of  the  unbiased  factor  for  the  biased 
c  approximate  MLE  of  the  scale  parameter  of  a  extreme  value  distribution 
c  with  a  sample  size  of  N. 


INTEGER  N,  S 

IF  (N.EQ.  10)  TK7M 
IF  (S.EQ.O) 

BIAS  =  0.  '•  39 
ELSEIF  (S.EQ.  i)  THEN 
BIAS  =  0.9275 
ELSEIF  (S.EQ.  2)  THEN 
BIAS  =  0. 9152 
ELSEIF  (S.EQ. 3)  THEN 
BIAS  =  0.9001 
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ELSEIF  CS.EQ.4)  THEN 
BIAS  =  0.8908 
ELSEIF  CS.EQ.8)  THEN 
BIAS  =  0. 8453 
ELSE 

BIAS  =  0. 7998 
ENDIF 

ELSEIF  (N.EQ. 15)  THEN 
IF  (S.EQ.O)  THEN 
BIAS  =  0. 95025 
ELSEIF  (S.EQ. 1)  THEN 
BIAS  =  0. 94715 
ELSEIF  (S.EQ. 2)  THEN 
BIAS  =  0. 94015 
ELSEIF  (S.EQ. 3)  THEN 
BIAS  =  0. 9309 
ELSEIF  (S.EQ  4)  THEN 
BIAS  =  0. y2495 
ELSEIF  (S.EQ. 8)  THEN 
BIAS  =  0. 89855 
ELSE 

BIAS  =  0.8718 
ENDIF 

ELSEIF  (N.EO. 20)  THEN 
IF  (S.EQ.O)  THEN 
BIAS  =  0. 9666 
ELSEIF  (S.EQ.  1)  THEN 
BIAS  =  0. 9668 
ELSEIF  (S.EQ.  2)  THEN 
BIAS  =  0.9651 
ELSEIF  (S.EQ.  3)  THEN 
BIAS  =  0. 9617 
ELSEIF  (S.EQ. 4)  THEN 
BIAS  =  0.9591 
ELSEIF  (S.EQ. 8)  THEN 
BIAS  =  C.9518 
ELSE 

BIAS  =  0. 9438 
ENDIF 
ENDIF 
RETURN 
END 


c  F.  Biased  APROXIMATE  MLE  Weibull  Shape  Parameter. 


SUBROUTINE  APROXMLE(Y,NN,A,BETAHAT) 

INTEGER  I,  NN,  R,  S,  A 

REAL*8  Yf*),  P(20),  Q(20),  ALPA(20),  BETA(20), 

+  BETAHAT,GAMKA,DELTA,B,C,D,E,M,BB,CC,DD,EE,MM.SIGMAHAT 


R  =  0 

S  =  NN  -  R  -  A 
DO  5  1=1,  A 
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P(I)  =  REAL(I)/(NN  +  1) 

Q(I)  =  1.0  -  P(I) 

ALPA(I)=1,0  +  L0G(Q(I))*(1.0  -  LOGC-LOG(Q(I)))) 

BETA(I)=  -LOG(Q(I)) 

5  CONTINUE 

GAMMA  =  -(Q(R+1)/P(R+1))*L0G(Q(R+1))*C1.0-L0GC-L0G(Q(R+1))))  + 

+  (QCR+1)/(P(R+1)**2))*((L0G(Q(R+1)))**2)*L0G(-L0G(Q(R+1))) 

DELTA  =  (Q(R+1)/P(R+1))*LOG(Q(R+1))*(1.0+(LOG(Q(R+1))/P(R+1))) 

BB  =  0.0 
CC  =  0.  0 
DD  =  0.  0 
EE  =  0.  0 
MM  =  0.0 

DO  10  I=R+1,  NN-S 
CC  =  CC  +  ALPA(I) 

BB  =  BB  +  BETA(I)*Y(I) 

MM  =  MM  +  BETACI) 

10  CONTINUE 

M  =  R*DELTA  +  S*BETA(NN-S)  +  MM 
B  =  (R*DELTA*y(R+l)  +  S*BETA(NN-S)-*YCNN"S)  +  BB)/M 
C  =  (R*GAMMA  -  S*(1,0  -  ALPA(NN-S))  +  CC)/M 

DO  15  I*=R+1,  NN-S 

DD  =  DD  +  ALPA(I)*(Y(I)  -  B) 

EE  =  EE  +  BETA(I)*((Y(I)  -  B)**2) 

15  CONTINUE 

D  =  R*GAMMA*CY(R+1)-B)  -  S*(l.  0-ALPA(NN-S))*(Y(NN-S)-B)  +  DD 
E  =  R*DELfiV*(CY(R+l)-B)**2)  +  S*BETA(NN-S)*((Y(NN-S)-B)**2)  +  EE 

SIGMAHAT  =(-D  +  SQRT(D’V*2  +  4.  0*E*FL0AT(A)))/(2.  0*FLOAT(A)) 
SIGMAHAT  =SIGMAHAT/BIAS(NN,S) 

BETAHAT  =  1.  0/SIGMAHAT 

RETURN 

END 


G.  Chi-Square  Quantile  Function. 


FUNCTION  CHISQD(P,N) 

Mod.-*zied  version  of  Algorithm  451  from  Communication  of  the  ACM 
August  1977  Vol.  16  No.  8. 

This  function  evaluates  the  quantile  at  the  probability  level  F 
(left  tail  area)  for  the  Chi-Square  Distribution  with 
N  degrees  of  freedom. 


REAL*8  P 
REAL  X 
INTEGER  IF 


DIMENSION  0(21),  AC  19) 


c 


DATA  0/  l,565326E-3 
+  1. 060438E-3 
+  -6. 950356E-3 
+  -1. 323293E-2 
+  2. 277679E-2 
+  -8. 986007E-3 
+  -1.513904E-2 
+  2.530010E-3 
+  -1.450117E-3 
+  5. 169654E-3 
+  -1. 153761E-2 
+  1. 128186E-2 
+  2. 607083E-2 
+  -0.2237368, 

+  9. 780499E-5 
+  -8.426S12E-4 
+  3. 125580E-3 
+  -8.553069E-3 
+  1. 348028E-4 
+  0.4713941, 

+  1.0000886  / 


i 


DATA  A/  1.264616E-2, 
+  -1.425296E-2, 

+  1.400483E-2, 

-5.886090E-3, 
T  -1. 091214E-2, 

+  -2.304527E-2, 

+  3. 135411E-3, 

+  -2.728484E-4, 

+  -9.699681E-3, 

+  1.316872E-2, 

+  2. 618914E-2, 

+  -0.2222222, 

+  5.406674E-5, 

+  3.483789E-5, 

+  -7. 274761E-4, 

+  3.292181E-3, 

+  -8. 729713E-3, 

+  0.4714045, 

+  1./ 


IF  (N-2)  10,  20,  30 
10  CALL  XFR0MP(.5*(1. -P),X.IF) 
CHISQD  =  X 

CHISQD  =  CHISQD*CHISQj 
RETURN 


c 

20  CHISQD  =  -2.*L0GC1. -P) 
RETURN 


30  F  =  N 

FI  =  1.  /F 

CALL  XFROMP(P,X,IF) 
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T  =  X 

F2  =  SQRT(F1)*T 

IF  (N.GE. (2+INT(4.*ABS(T))))  GO  TO  40 
c 

CHISQD  =  (((((( (C(1)*F2+C(2))*F2+C(3))*F2+C(4))*F2 
+  +C( 5 ) )*F2+C( 6) )*F2+C( 7 ) )*Fl+( ( ( ( ( ( C( 8)+C( 9 )*F2)*F2 

+  +C( 10) )*F2+C( 11) )*F2+C( 12) )*F2+C( 13) )*F2+C( 14) ) )*F1+ 

+  C ( ( ( (C( 15)*F2+C( 16) )*F2+C( 17 ) )*F2+C( 18) )*F2 

+  +C(19))*F2+CC20))*F2+C(21) 

c 

GO  TO  SO 
c 

40  CHISQD  =  (((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 

+  +A(5) )*F2+A( 6) ) )*Fl+( ( ( ( ( A( 7)+A( 8)*F2)*F2+A( 9) )*F2 

+  +A( 10) )*F2+A( 11) )*F2+A( 12) ) )*Fl+( ( ( ( ( A( 13)*F2 

+  +A( 14) )*F2+A( 15) )*F2+A( 16) )*F2+A( 17) )*F2*F2 

+  +A(18))*F2+A(19) 

c 

50  CHISQD  =  CHISQD*CHISQD*CHISQD*F 
C 

RETURN 

END 


c  H.  Standard  Normal  Variate  Computation  Subroutine. 


SUBROUTINE  XFROMPC  P , X , IFAULT) 
c 

c  Algorithm  AS  24  J.  R.  STAT. SOC.  C  (1969)  Vol.  18.  No.  3. 
c 

c  This  subroutine  compute  the  standard  normal  deviate  X  for  the 
c  specified  left  tail  area  P. 


REAL*8  P 
DIMENSION  A(5) 

DIMENSION  CONNOR  (17),  HSTNGS(6) 


DATA  CONNOR/  8. 0327350124E-17, 
1.4483264644E-15, 
2.4668270103E-14, 
3. 9554295 164E- 13, 
5.9477940136E-12, 
8. 3507027951E-11, 
1.0892221037E-9, 
1.3122532964E-8, 
1.4503852223E-7, 
1.4589169001E-6, 
1. 3227513228E-5, 
1. 0683760684E-4, 
7.5757575758E-4, 
4. 6296296296E-3, 
2.3809523810E-2, 
0.1, 

0.33333333333  / 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


c 
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c 

c 

c 


c 


c 


c 

c 

c 


c 


c 


FATA  RTHFPI  /  1. 2S331A1373  / 


DATA  RRT2PI  /  0. 3989A22804  / 
DATA  TERMIN  /  l.OE-11  / 


DATA  HSTNGS 
+ 

+ 

+ 

+ 

+ 


2.515517, 
0. 802853, 
0.010328, 
1. A32788, 
0. 189269, 
0.001308  / 


IFAULT  =  1 

IF  CCP.LE.  0.0).OR.  (P.6E.  1.0))  GO  TO  100 
IFAULT  =  0 


Get  first  approximation  Xo  to  deviate  by  Hasting’s  formula. 
B  =  P 

IF  (B.GT.  0.5)  B  =  1.0  "  B 


F  =  -  LOG(B) 

E  =  SQRT(F+F) 

XO  =  -E  +  ((HSTNGS(3)*E+HSTNGS(2))*E+HSTNGS(1))/ 

+  ((CHSTNGS(6)*E+HSTNGS(5))*E+HSTNGS(A))*E+1. 0) 


IF  (XO.  LT.  0. 0)  GO  TO  1 
XO  =  0.  0 
PO  =  0.  5 
XI  =  -RTHFPI 
GO  TO  7 


c  Find  the  area  PO  corresponding  to  XO 
c 

1  Y  =  X0**2 

IF  (XO.LE.  -1.9)  GO  TO  3 
Y  =  -0.  5*Y 
c 

c  (1)  series  approximation 
c 

PO  =  CONNOR(l) 

DO  2  L  -  2,  17 

2  PO  =  P0*Y  +  CONNOR(L) 

PO  =  (P0*y+1.  0)*X0 

XI  =  -(PO+RTHFPI)*EXP(-Y) 

PO  =  PO*RRT2PI  +0.5 
GO  TO  7 
c 

c  (2)  continued  fraction  approximation 
c 

3  Z  =  1.  O/Y 
A(2)  =  1.0 
A(3)  =  1.0 
A(A)  =  Z  +  1.0 
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A(5)  =  1.0 
W  =  2.0 
c 

4  DO  6  L  =  1,  3,  2 

DO  5  J  =  1,  2 
K  =  L  +  J 
KA=  7  -  K 

5  A(K)  =  ACKA)  +  A(K)*W*Z 

6  W  =  W  +  1.  0 
c 

APPRXU  =  A(2)/A(3) 

APPRXL  =  A(5)/A(4) 

C  =  APPRXU  -  APPRXL 
c 

IF  CC.GE.TERMIN)  GO  TO  4 
c 

XI  =  APPRXL/XO 
PO  =  -X1*RRT2PI*EXPC-0.5*Y) 
c 

c  Get  accurate  value  of  deviate  by  Taylor  Series, 
c  (XI,  X2,  X3  are  derivatives  for  the  Taylor  series) 
c 

7  D  =  F  +  LOG(PO) 

X2  =  X0*X1*X1  -  XI 
X3  =  Xl**3  +  2. 0*X0*X1*X2  -  X2 
X  =  ((X3*D/3.  0+X2)*D/2. 0+Xl)*D  +  XO 
c 

IF  (P.LE.  0.5)  GO  TO  100 
X  =  -X 
100  RETURN 
END 
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APPENDIX  D. 


TABULATED  RUN  RESULTS  FOR  RETP 
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Table  lA  :  8  Exp  in  Series,  RS  =  0.9305  (Hi) 


min  ?.  =  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

- i 

Measures  of  accuracy 

Plan 

Freedom 

cc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.9262 

0.9630 

10  failed 

(160) 

0.2 

0.9256 

0.9220 

NFC  =80 

2(NFC+ 

0.1 

0.9196 

0.9970 

NCOMP) 

(176) 

0.2 

0.9188 

0.9910 

2NFC- 

0.1 

!  0.9296 

0.9209 

NCO.MP 

(152) 

0.2 

0.9291 

0.8430 

2(NFC- 

0.1 

0.9329 

0.8400 

NCOMP) 

(144) 

0.2 

0.9325 

0.7300 

2 

Test  15  until 

2NFC 

0.1 

0.9277 

0.9550 

15  failed 

(240) 

0.2 

0.9273 

0.9080 

NFC  =120 

2(NFC  + 

0.1 

0.9233 

0.9899 

NCOMP) 

(256) 

0.2 

0.9228 

0.9750 

2NFC- 

0.1 

0.9298 

0.9159 

NCOMP 

(232) 

0.2 

0.9296 

0.8329 

2(NFC- 

0.1 

0.9321 

0.8440 

NCOMP) 

(224) 

0.2 

0.9318 

0.7470 

3 

Test  15  until 

2NFC 

0.1 

0.9268 

0.9550 

1 1  failed 

(176) 

0.2 

0.9262 

0.9159 

NFC  =88 

2(NFC+ 

0.1 

0.9208 

0.9960 

NCOMP) 

(192) 

0.2 

0.9200 

0.9880 

2NFC- 

0.1 

0.9298 

0.9159 

NCOMP 

(168) 

0.2 

0.9293 

0.8430 

2(NFC- 

0.1 

0.9228 

0.8430 

NCOMP) 

(160) 

0.2 

0.9324 

0.7350 

Table  lA  :  8  Exp  in  Series,  RS  =  0.9305  (Hi)  (Cont...) 

_ min  A  =  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC =56 

2NFC 

(112) 

0.1 

0.9241 

0.9700 

0.2 

0.9229 

0.9310 

2(NFC+ 
NCOMP) 
(128)  ■ 

0.1 

0.9145 

0  9980 

0.2 

0.9130 

0.9940 

2NFC- 

NCOMP 

(104) 

0.1 

0.9289 

0.9190 

0.2 

0.9280 

0.8530 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.9338 

0.8350 

0.2 

0.9331 

0.7200 

5 

Test  15  until 

3  failed 

NFC =24 

2NFC 

(48) 

0.1 

0.9145 

0.9859 

0.2 

0.9129 

0.9750 

2{NFC+ 

NCOMF 

(64) 

0.1 

0.8907 

1.0000 

0.2 

0.8877 

1.0000 

2NFC- 

NCOMP 

(40) 

0.1 

0.9268 

0.9439 

02 

0.9260 

0.8600 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.9394 

0.7530 

0.2 

0.9394 

0.6339 
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Table  IB  :  8  Exp  in  Series,  RS  =  0.8025  (Low) 


min  =  0.0010  f/hr,  max  A  =  0.0045  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

oc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.9710 

10  failed 

(160) 

0.2 

0.9310 

NFC  =80 

2(NFC  + 

0.1 

0.7735 

0.9990 

NCOMP) 

tl76) 

0.2 

0.7707 

0.9920 

2n:'.'- 

0.1 

0.7994 

0.9299 

NCOMP 

(152) 

0.2 

0.7975 

0.8500 

2(NFC- 

0.1 

0.8083 

0.8480 

NCOMP) 

(144) 

0.2 

0.8066 

0.7410 

2 

Test  15  until 

2NFC 

0.1 

0.7940 

0.9620 

15  failed 

(240) 

0.2 

0.7933 

0.9159 

NFC  =120 

2(NFC+ 

0.1 

0.7826 

0.9930 

NCOMP) 

(256) 

0.2 

0.7816 

0.9809 

2NFC- 

0.1 

0.'.78 

0.9240 

NCOMP 

(232) 

0.2 

0./992 

0.8400 

2(NFC- 

0.1 

0.8057 

0.8580 

NCO.MP) 

(224) 

0.2 

0.8052 

0.7550 

3  1 

Test  15  until 

2NFC 

0.1 

0.7919 

0.9660 

1 

11  failed 

(176)  j 

0.2 

0.7896 

0.9209 

i 

NFC  =88 

2(NFC+  i 

0.1 

0.7764 

0.9960 

i 

1 

NCOMP)  1 
(192) 

0.2 

0.7735 

0.9890 

2NFC- 

0.1 

0.7999 

0.9220 

1 

NCOMP 

(168) 

0.2 

0.7978 

0.8550 

2(NFC- 

0.1 

0.8079 

0.8550 

NCOMP) 

(160) 

0.2 

0.8061 

0.746C 
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Table  IB  :  8  Exp  in  Series,  RS  =  0.8025  (Low)  (Cent...) 

_ min  /.  =  0.0002  f/hr,  max  2  =  0.0016  f/hr.  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC  =  56 

2NFC 

(112) 

0.1 

0.7851 

0.9740 

0.2 

0.7816 

0.9380 

2(NFC+ 

NCOMP) 

(128) 

0.1 

0.7605 

0.9980 

0.2 

0.7559 

0.9960 

2NFC- 

NCOMP 

(104) 

0.1 

0.7978 

0.9250 

0.2 

0.7948 

0.8600 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.8107 

0.8410 

0.2 

0.8083 

0.7270 

5 

Test  15  until 

3  failed 

NFC =24 

2NFC 

(48) 

0.1 

0.7593 

0.9890 

0.2 

0.7555 

0.9770 

2(NFC  + 
NCOMP) 
(64) 

0.1 

0.7001 

1.0000 

0.2 

0.6928 

1.0000 

2NFC- 

NCOMP 

(40) 

0.1 

0.7912 

0.9489 

0.2 

0.7892 

0.8720 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.8249 

0.7630 

0.2 

0.8248 

0.6420 
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Table  2A  :  8  Wei  in  Series,  RS  =  0.9798 


Deg.  of 
Freedom 

a 

2NFC 

0.1 

(160) 

0.2 

2(NFC  + 

0.1 

NCOMP) 

(176) 

0.2 

2NFC- 

0.1 

NCOMP 

(152) 

0.2 

2(NFC- 

0.1 

NCOMP) 

(144) 

0.2 

2NFC 

0.1 

(240) 

0.2 

2(NFC 

0.1 

NCOi\.:n 

(256) 

0.2 

2NFC- 

0.1 

NCOMP 

(232) 

0.2  ' 

2(NFC- 

0.1 

NCOMP) 

(224) 

0.2 

2NFC 

0.1 

(176) 

0.2 

2(NFC  + 

0.1 

jCOMP) 

(192) 

0.2 

2NFC- 

0.1 

NCOMP 

(168) 

0.2 

2(NFC- 

0.1 

NCOMP) 

(160) 

0.2 

Measures  of  accuracy 


RSLOW 


0.9602 


0.9506 


0.9565 


0.946i 


C.9619 


0.9529 


0.9638 


0.9553 


0  9706 


0.9649 


0.9687 


0.9627 


0.9715 


0.9659 


0.9724 


0.9671 


0.9  "79 


0.9606 


0.9716 


0.9652 


0.9729 


0.9668 


LEVEL 


0.9859 


0.9820 


0.9890 


0.9870 


0.9820 


0.9809 


0.9800 


0.9770 


0.9770 


0.9719 


0.9820 


0.9790 


0.9730 


0.9640 


0.9650 


0.9590 


0.9719 


0.9650 


0.9800 


0.9760 


0.9660 


0.9590 


0.9590 


0.9510 


Table  2A  :  8  Wei  in  Series,  RS  =  0.9798  (Hi)  (Cent...) 

min  A  --  0.001  f/hr,  max  X  -  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC  =  56 

2NFC 

(112) 

0.1 

0.9764 

0.9400 

0.2 

0.9668 

0.9280 

2(NFC  + 
NCOMP) 
(128) 

0.1 

0.9733 

0.9579 

0.2 

0.9624 

0.9510 

2NFC- 

NCOMP 

(104) 

0.1 

0.9779 

0.9240 

0.2 

0.9691 

0.9110 

2(NrC- 

NCOMP) 

(96) 

0.1 

0.9795 

0.9080 

0.2 

0.9713 

0.8900 

5 

Test  15  until 

3  failed 

NFC  =  24 

2NFC 

(48) 

0.1 

0.9889 

0.8000 

0.2 

0.9814 

0.7770 

2(NFC+ 

NCO.MP) 

(64) 

0.1 

0.9857 

0.8469 

0.2 

0.9758 

0.8360 

2NFC- 

NCOMP 

(40) 

0.1 

0.9906 

0.7510 

0.2 

0.9843 

0.7280 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.9922 

0.6890 

0.2 

0.9872 

0.6530 

Table  2B  :  8  Wei  in  Series,  RS  =  0.8323  (Low) 

_ min  X  =  0.003  f/hr,  max  2  =  0.024  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

2NFC 

(160) 

0.1 

0.7616 

0.9800 

0.2 

0.7383 

0.9740 

2(NFC+ 

NCOMP) 

(176) 

0.1 

0.7425 

0.9880 

0.2 

0.7171 

0.9840 

2NFC- 

NCOMP 

(152) 

0.1 

0.7714 

0.9730 

0.2 

0.7491 

0.9680 

2(NFC- 

NCOMP) 

(140 

0.1 

0.7813 

O.9680 

0.2 

0.7601 

0.9620 

2 

Test  15  until 

15  failed 

NFC=  120 

2NFC 

'10) 

0.1 

0.8002 

0.9640 

0.2 

0.7766 

0.9529 

2(NFC+ 

NCO.MP) 

(256) 

0.1 

0.7891 

0.9730 

0.2 

0.7641 

0.9669 

2NFC- 

NCOMP 

(232) 

0.1 

0.9059 

0.9560 

0.2 

0.7829 

0.9430 

2(NFC- 

NCOMP) 

(224) 

0.1 

0.8116 

0.9450 

0.2 

0.7893 

0.9320 

3 

Test  15  until 

11  failed 

NFC  =88 

2NFC 

(176) 

0.1 

0.8056 

0.9529 

0.2 

0.7803 

0.9430 

2(NFC+ 

NCOMP) 

(192) 

0.1 

0.7909 

0.9710 

0.2 

0.7636 

0.9640 

2NFC- 

NCOMP 

(168) 

0.1 

0.8131 

0.9430 

0.2 

0.7888 

0.9270 

2(NFC- 

NCOMP) 

(160) 

0.1 

0.8206 

0.9270 

0.2 

0.7994 

0.9060 
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nrs. 


Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

(X 

RSLOW 

LEVEL 

Test  15  until 

2NFC 

0.1 

0.8296 

0.9030 

7  failed 

(112) 

0.2 

0.8064 

0.8830 

NFC =56 

2(NFC+ 

0.1 

0.8095 

0.9470 

NCOMP) 

(128) 

0.2 

0.7832 

0.9320 

2NFC- 

0.1 

0.8399 

0.8800 

NCOMP 

(104) 

0.2 

0.8183 

0.8540 

2(NFC- 

0.1 

0.8504 

0.8480 

NCOMP) 

(96) 

0.2 

0.8304 

0.8070 

Test  15  until 

2NFC 

0.1 

0.8801 

0.7449 

3  failed 

(48) 

0.2 

0.8569 

0.7060 

NFC =24 

2(NFC+ 

0.1 

0.8475 

0.8660 

NCOMP) 

(64) 

0.2 

0.8171 

0.8410 

2NFC. 

0.1 

0.8970 

0.6540 

NCOMP 

(40) 

0.2 

0.8778 

0.5969 

2(NFC- 

0.1 

0.9145 

0.5270 

NCOMP) 

(32) 

0.2 

0.8994 

0.4470 

Table  3 A  :  4  Exp  and  4  Wei  (mixed)  in  Series,  RS  =  0.9801  (Hi) 
_ min  /  =  0.0002  f/hr,  max  X  =  0.0008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

Measures  of  accuracy 

Plan 

Freedom 

oc 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.9796 

0.9190 

10  failed 

(160) 

0.2 

0.9791 

0.8720 

NFC  =80 

2(NFC+ 

0.1 

0.9777 

0.9740 

NCOMP) 

(176) 

0.2 

0.9771 

0.9489 

2NFC- 

0.1 

0.9805 

0.8720 

NCOMP 

(152) 

0.2 

0.9801 

0.7980 

2(NFC- 

0.1 

0.9815 

•  0.7950 

NCOMP) 

(144) 

0.2 

0.9811 

0.7000 

2 

Test  15  until 

2NFC 

0.1 

0.9802 

0.8950 

15  failed 

(240) 

0.2 

0.9799 

0.8250 

NFC  =120 

2(NFC+ 

0.1 

0.9789 

0.9510 

NCOMP) 

(256) 

0.2 

0.9786 

0.9209 

2NFC- 

0.1 

0.9808 

0.8370 

NCOMP 

(232) 

0.2 

0.9805 

0.7530 

2(NFC- 

0.1 

0.9814 

0.7670 

NCOMP) 

(224) 

0.2 

0.9812 

0.6670 

3 

Test  15  until 

2NFC 

0.1 

0.9802 

0.8920 

11  failed 

(176) 

0.2 

0.9797 

0.8400 

NFC  =88 

2(NFC+ 

0.1 

0.9785 

0.9590 

NCOMP) 

(192) 

0.2 

0.9779 

0.9320 

2NFC- 

0.1 

0.9810 

0.8410 

NCOMP 

(168) 

0.2 

0.9806 

0.7550 

2(NFC- 

0.1 

0.9819 

0.7580 

NCOMP) 

(160) 

0.2 

0.9815 

0.6590 

Table  3A  :  4  Exp  and  4  Wei  in  Series,  RS  =  0.9801  (Hi)  (Cont...) 
_ min  ).  =  0.0002  f/hr,  max  X  =  0.0008  f/hf,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC =56 

2NFC 

(112) 

0.1 

0.9803 

0.8950 

0.2 

0.9797 

0.8260 

2(NFC+ 

NCOMP) 

(128) 

0.1 

0.9777 

0.9690 

0.2 

0.9770 

0.9439 

2NFC- 

NCOMP 

(104) 

0.1 

0.9815 

0.8060 

0.2 

0.9811 

0.7300 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.9828 

0.7090 

0.2 

0.9825 

0.6000 

5 

Test  15  until 

3  failed 

NFC =24 

2NFC 

(48) 

0.1 

0.9815 

0.8730 

0.2 

0.9800 

0.8020 

2(NFC+ 

NCOMP) 

(64) 

0.1 

0.9757 

0.9740 

0.2 

0.9739 

0.9560 

2NFC- 

NCOMP 

(40) 

0.1 

0.9839 

0.7250 

0.2 

0.9831 

0.6260 

2(NFC- 

NCOiMP) 

(32) 

0.1 

0.9868 

0.5150 

0.2 

0.9862 

0.4140 
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Table  3B  :  4  Exp  and  4  Wei  (mixed)  in  Series,  RS  =  0.8089  (Low) 


min  A  =  0 

.002  f/hr,  max  A  = 

0.008  f/hr,  UT 

=  10  hrs. 

S/N 

Test 

Deg.  of 

a 

.Measures  of  accuracy 

Plan 

Freedom 

RSLOW 

LEVEL 

1 

Test  10  until 

2NFC 

0.1 

0.7957 

0.9570 

10  failed 

(160) 

0.2 

0.7914 

0.9240 

NFC  =80 

2(NFC+ 

0.1 

0.7789 

0.9920 

NCOMP) 

(176) 

0.2 

0.7738 

0.9770 

2NFC- 

0.1 

0.8042 

0.9240 

NCOMP 

(152) 

0.2 

0.8003 

0.8730 

2(NFC- 

0.1 

0.8129 

0.8710 

NCOMP) 

(144) 

0.2 

0.8093 

0.7950 

2 

Test  15  until 

2NFC 

0.1 

0.8033 

0.9340 

15  failed 

(240) 

0.2 

0.8007 

0.8929 

NFC=  120 

2(NFC+ 

0.1 

0.7923 

0.9740 

NCOMP) 

(256) 

0.2 

0.7893 

0.9500 

2NFC- 

0.1 

0.8089 

0.9010 

NCOMP 

(232) 

0.2 

0.8064 

0.8260 

2(NFC- 

0.1 

0.8145 

0.8400 

NCOMP) 

(224) 

0.2 

0.8122 

0.7750 

3 

Test  15  until 

2NFC 

0.1 

0.8045 

0.9290 

11  failed 

(176) 

0.2 

0.7993 

0.8810 

NFC  =88 

2(NFC+ 

0.1 

0.7897 

0.9850 

NCOMP) 

(192) 

0.2 

0.7839 

0.9570 

2NFC- 

0.1 

0.8119 

0.8820 

NCOMP 

(168) 

0.2 

0.8072 

0.8189 

2(NFC- 

0.1 

0.8195 

0.8200 

NCOMP) 

(160) 

0.2 

0.8151 

0.7300 
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Table  3B  ;  4  Exp  and  4  Wei  in  Series,  RS  =  0.8089  (Low)  (Cont...) 

rrvin  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed 

NFC =56 

2NFC 

(112) 

0.1 

0.8056 

0.9140 

0.2 

0.8017 

0.8600 

2(NFC+ 

NCOMP) 

(128) 

0.1 

0.7830 

0.9760 

0.2 

0.7780 

0.9610 

2NFC- 

NCOMP 

(104) 

0.1 

0.8172 

0.8460 

0.2 

0.8139 

0.7630 

2(NFC- 

NCOMP) 

(96) 

0.1 

0.8290 

0.7410 

0.2 

0.8262 

0.6490 

5 

Test  15  until 

3  faUed 

NFC =24 

2NFC 

(48) 

0.1 

0.8174 

0.8670 

0.2 

0.8085 

0.8060 

2(NFC+ 

NCOMP) 

(64) 

O.l 

0.7703 

0.9780 

0.2 

0.7571 

0.9550 

2NFC- 

NCOMP 

(40) 

0.1 

0.8424 

0.7219 

0.2 

0.8357 

0.6050 

2(NFC- 

NCOMP) 

(32) 

0.1 

0.8685 

0.4859 

0.2 

0.8641 

0.3810 
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Table  4A  :  8  Exponential  in  Parallel,  RS  =  0.9345  (Hi) 

min  X  =  0.1000  f/hr,  max  X  =  0.1525  fyhr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

160 

0.1 

0.9306 

0.9290 

0.2 

0.9320 

0.8180 

2 

Test  15  until 

15  failed 

NFC  =120 

240 

0.1 

0.9325 

0.9190 

0.2 

0.9324 

0.8210 

3 

/ 

Test  15  until 

11  failed 

NFC  =88 

176 

0.1 

0.9309 

0.9240 

0.2 

0.9312 

0.8329 

4 

Test  15  until 

7  failed 

NFC  =  56 

112 

0.1 

0.9307 

0.9209 

0.2 

0.9308 

0.8310 

5 

Test  15  until 

3  failed 

NFC =24 

48 

0.1 

0.9280 

0.9220 

0.2 

0.9286 

0.8370 

1 
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Table  4B  ;  8  Exponential  in  Parallel,  RS  = 
min  X  =  0.1000  f/hr,  max  /  = 


0.8262  (Low) 

0.2575  f/hr,  UT  =  10  hrs. 


S/N 

Test 

Deg.  of 

a 

Measures  of  accuracy 

Plf-n 

Freedom 

RSLOW 

LEVEL 

1 

1  est  i'  aiitil 

160 

0.1 

0.8309 

0.8870 

10  m;-.ed 

NFC- 80 

0.2 

0.8^.,8 

0.7880 

2 

Test  15  V  'til 

240 

0.1 

0.8290 

0.8870 

15  failed 

NFC- 120 

0.2 

0.8283 

0.78C0 

1 

Test  15  until 

176 

0.1 

0.8266 

0.8990 

U  failed 

F  :  ■-  88 

0.2 

0.8264 

0.7990 

j  4 

■>  a  a.  15  until 

112 

0.1 

0.8300 

0.8860 

7  failed 

NFC- 56 

0.2 

0.8284 

0.7890 

Test  15  until 

48 

0.1 

0.8354 

0.8810 

3  failed 

NFC =24 

0.2 

0.8364 

0.7700 

Table  5A  :  8  Wei  in  Parallel  ,  RS  ==  0.9265  (Hi) 

min  X  =  0.100  f/hr,  max  X  =  0.128  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  Q.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =  80 

161 

0.1 

0.9154 

0.9540 

0.2 

0.9157 

0.8880 

2 

Test  15  until 

15  failed 

NFC=  120 

241 

O.I 

0.9179 

0.9510 

0.2 

0.9170 

0.8770 

3 

Test  15  until 

11  failed 

NFC  =  88 

177 

0.1 

0.9248 

0.9180 

0.2 

0.9250 

0.8130 

4 

Test  15  until 

7  failed 

NFC  =  56 

113 

0.1 

0.9323 

0.8630 

0.2 

0.9318 

0.7589 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.9152 

0.9270 

0.2 

0.8980 

0.8850 
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Table  5B  :  8  Wei  in  Parallel  ,  RS  =  0.8351  (Low) 

min  A  =0.100  f/hr,  max  A  =  0.163  f)'hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

N1'C=80 

161 

0.1 

0.8296 

0.9180 

0.2 

0.8331 

0.8140 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.8322 

0.9100 

0.2 

0.8314 

0.8279 

3 

Test  15  until 

11  failed 

NP'C=88 

177 

0.1 

0.8494 

0.8380 

0.2 

0.8489 

0.7130 

4 

Test  15  until 

7  failed 

NFC =56 

113 

0.1 

0.8713 

0.7560 

0.2 

0.8705 

0.6450 

5 

Test  15  until 

3  failed 

NFC  =24 

51 

0.1 

0.8678 

0.8360 

0.2 

0.8469 

0.7770 

10 


Table  6A  :  4  EXP  and  4  Wei  (mixed)  in  Parallel,  RS  =  0.9408  (Hi) 
min  X  -  0.100  f/hr,  max  X  =  0.130  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

.Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.9342 

0.9420 

0.2 

0.9366 

0.8500 

2 

Test  15  until 

15  failed 

NFC=  120 

241 

0.1 

0.9362 

0.9330 

0.2 

0.9363 

0.8510 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9389 

0.9240 

0.2 

0.9378 

0.8310 

4 

Test  15  until 

7  failed 

NFC =56 

113 

0.1 

0.9416 

0.8920 

0.2 

0.9427 

0.7819 

5 

Test  15  until 

3  failed 

NFC  =  24 

51 

0.1 

0.9383 

0.9080 

0.2 

0.9325 

0.8410 

\ 
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Table  6B  :  4  EXP  and  4  Wei  (mixed)  in  Parallel,  RS  =  0.8170  (Low) 
min  A  =  0.100  f/hr,  max  A  =  0.220  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

.Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =  80 

161 

0.1 

0.8264 

0.8550 

0.2 

0.8273 

0.7460 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.8262 

0.8630 

0.2 

0.8236 

0.7560 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.8333 

0.8290 

0.2 

0.8334 

0.7150 

4 

Test  15  until 

7  failed 

NFC =56 

113 

0.1 

0.8488 

0.7660 

0.2 

0.8469 

0.6470 

5 

Test  15  until 

3  failed 

NFC =24 

51 

0.1 

0.8650 

0.7819 

0.2 

0.8563 

0.6709 
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Table  7A  :  8  Exp  in  Series-  Parallel,  RS  =  0.9249  (Hi) 

min  X  =  0.0003  f/hr,  max  X  =  0.0024  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =2NFC _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

.Measures  of  accura':y 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

160 

0.1 

0.9241 

0.9159 

0.2 

0.9226 

0.8620 

2 

Test  15  until 

15  failed 

NFC  =120 

240 

0.1 

0.9253 

0.8960 

0.2 

0.9236 

0.8360 

3 

Test  15  until 

11  failed 

NFC  =  88 

176 

0.1 

0.9252 

0.8950 

0.2 

0.9226 

0.8390 

4 

Test  15  until 

7  failed 

NFC  =  56 

112 

0.1 

0.9222 

0.9260 

0.2 

0.9204 

0.8789 

5 

Test  15  until 

3  failed 

NFC =24 

48 

0.1 

0.9153 

0.9620 

0.2 

0.9123 

0.9260 
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Table  7B  :  8  Exp  in  Series-  Parallel,  RS  =  0.8228  (Low) 

min  X  =0.00075  f/hr,  max  X  =0.0060  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC _ _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

I 

Test  10  until 

10  failed 

NFC  =  80 

160 

0.1 

0.8210 

0.9159 

0.2 

0.8175 

0.8620 

2 

Test  15  until 

15  failed 

NFC  =120 

240 

0.1 

0.8237 

0.8960 

0.2 

0.8199 

0.8360 

3 

Test  15  until 

11  failed 

NFC  =88 

176 

0.1 

0.8234 

0.8950 

0.2 

0.8176 

0.8390 

4 

Test  15  until 

7  failed 

NFC =56 

112 

0.1 

0.8168 

0.9260 

0.2 

0.8127 

0.8789 

5 

Test  15  until 

3  failed 

NFC  =  24 

48 

0.1 

0.8015 

0.9620 

0.2 

0.7949 

0.9260 
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Table  8A  :  8  Wei  in  Series-  Parallel,  RS  =  0.9328  (Hi) 

min  X  =  0.002  f/hr,  max  /  =  0.016  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.9106 

0.9579 

0.2 

0.8939 

0.9510 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.9247 

0.9330 

0.2 

0.9108 

0.9230 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9266 

0.9220 

0.2 

0.9149 

0.9119 

4 

Test  15  until 

7  failed 

NFC  =56 

113 

0.1 

0.9414 

0.8609 

0.2 

0.9250 

0.8419 

5 

Test  15  until 

3  failed 

NFC  =24 

51 

0.1 

0.9699 

0.7079 

0.2 

0.9536 

0.6740 
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Table  8B  :  8  Wei  in  Series-  Parallel,  RS  =  0.8321  (Low) 

min  X  —  0.00325  f/hr,  max  X  =  0.026  f/hr,  UT  =  10  hrs. 


Deg.  of 
Freedom 


a 

Measures  of  accuracy 

RSLOW 

LEVEL 

0.1 

0.7963 

0.9560 

0.2 

o.im 

0.9470 

0.1 

0.8204 

0.9310 

0.2 

0.7986 

0.9200 

0.1 

0.8250 

0.9190 

0.2 

0.8046 

0.9000 

0.1 

0.8485 

0.8590 

0.2 

0.8237 

0.8329 

0.1 

0.8999 

0.7500 

0.2 

0.8599 

0.6950 

16 
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Table  9A  :  4  EXP  and  4  Wei  in  Series-  Parallel,  RS  =  0.9276  (Hi) 
min  X  =  0.005  f/hr,  max  X  =  0,020  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =  80 

161 

0.1 

0.9096 

0.9400 

0.2 

0.8904 

0.9349 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.9238 

0.9170 

0.2 

0.9079 

0.9040 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.9266 

0.9050 

0.2 

0.9109 

0.8850 

4 

Test  15  until 

7  failed 

NFC  =  56 

113 

0.1 

0.9393 

0.8550 

0.2 

0.9203 

0.8390 

5 

Test  15  until 

3  failed 

NFC =24 

51 

0.1 

0.9637 

0.7430 

0.2 

0.9451 

0.7200 

n 


'1 
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Table  9B  :  4  EXP  and  4  Wei  in  Series-  Parallel,  RS  =  0.8248  (Low) 
min  X  =  0.008  f/hr  max  X  =  0.032  f/hr,  UT  =  10  hrs. 
_ Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) _ 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

OL 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =80 

161 

0.1 

0.7909 

0.9420 

0.2 

0.7627 

0.9360 

2 

Test  15  until 

15  failed 

NFC  =120 

241 

0.1 

0.8164 

0.9150 

0.2 

0.7915 

0.9020 

3 

Test  15  until 

11  failed 

NFC  =88 

177 

0.1 

0.8236 

0.9010 

0.2 

0.7964 

0.8820 

4 

Test  15  until 

7  failed 

NFC =56 

113 

0.1 

0.8412 

0.8660 

0.2 

0.8119 

0.8419 

5 

Test  15  until 

3  failed 

NFC  =24 

51 

0.1 

0.8683 

0.8170 

0.2 

0.8289 

0.7940 
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Table  lOA  :10  EXP  and  5  Wei  in  Series-  Parallel,  RS  =  0.9472  (Hi) 
Exp  :  min  X  =  0.025  f/hr,  max  X  =  0.075  f/hr. 

Wei  :  min  X  =  0.055  f/hr,  max  X  =  0.095  f/hr. 

UT  =  10  hrs. 


Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

a 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =150 

301 

0.1 

0.9383 

0.9710 

0.2 

0.9369 

0.9550 

2 

Test  20  until 

20  failed 

NFC  =300 

601 

0.1 

0.9449 

0.9370 

0.2 

0.9432 

0.8950 

3 

Test  20  until 

15  failed 

NFC  =225 

451 

0.1 

0.9444 

0.9349 

0.2 

0.9430 

0.8979 

4 

Test  20  until 

10  failed 

NFC  =150 

301 

0.1 

0.9435 

0.9290 

0.2 

0.9415 

0.9000 

5 

Test  20  until 

5  failed 

NFC =75 

152 

0.1 

0.9344 

0.9809 

0.2 

0.9310 

0.9520 

1 
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Table  lOB  :10  EXP  and  5  Wei  in  Series-  Parallel,  RS  =  0.8324  (Low) 
Exp  :  min  /  =  0.025  f/hr,  max  X  =  0.175  f/hr. 

Wei  :  min  X  =  0.055  f/hr,  max  X  =  0.125  f/hr. 

UT  =  10  hrs. 


Degrees  of  Freedom  =  2NFC  +  0.5(NC/NF) 


S/N 

Test 

Plan 

Deg.  of 
Freedom 

. 

a. 

Measures  of  accuracy 

RSLOW 

LEVEL 

1 

Test  10  until 

10  failed 

NFC  =150 

301 

0.1 

0.8103 

0.9669 

0.2 

0.8081 

0.9430 

2 

Test  20  until 

20  failed 

NFC  =  300 

601 

0.1 

0.8247 

0.9380 

0.2 

0.8219 

0.8990 

3 

Test  20  until 

15  failed 

NFC =225 

451 

0.1 

0.8236 

0.9470 

0.2 

0.8204 

0.8990 

4 

Test  20  until 

10  failed 

NFC  =150 

301 

0.1 

0.8232 

0.9500 

0.2 

0.8166 

0.8950 

5 

Test  20  until 

5  failed 

NFC =75 

152 

0.1 

0.8008 

0.9760 

0.2 

0.7917 

0.9489 

0 


APPENDIX  E.  TABLE  OF  CHI-SQUARE  DISTRIBUTION 

Table  21  in  this  appendix  provides  the  eightieth  and  ninetieth 
percentile  points  for  the  chi-square  distribution  with  degrees  of  freedom 
ranging  from  1  to  402  in  increments  of  one.  The  percentile  points  appear 
under  the  column  headed  CHI  in  Table  21. 

The  computational  algorithm  used  to  construct  these  percentile  points 
is  defined  in  the  Chi-Square  Quantile  Function  routine  together  with  the 
Standard  Normal  Variate  Computation  routine  on  pages  79  -  83  in  this 
thesis. 

Values  of  these  chi-square  percentile  points  for  degrees  of  freedom  as 
large  as  601  were  needed  in  this  thesis.  The  extensive  set  of  chi-square 
percentile  values  in  Table  21  is  provided  as  a  convenience  for  a  reader 
who  may  want  to  apply  this  methodology  to  a  particular  system. 
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Table  21  :  Chi-square  distribution 


C" - - - - - - - - - - 

Left  tail  area  =  0.  9  Left  tail  area  =  0.  8 


c 


DF 

CHI 

DF 

CHI 

DF 

CHI 

DF 

CHI 

DF 

CHI 

DF 

CHI 

1 

2.706 

47 

59. 774 

93 

110. 850 

1 

1.642 

47 

54. 905 

93 

104. 241 

2 

4.  605 

48 

60.  907 

94 

111. 944 

2 

3.219 

48 

55. 993 

94 

105. 303 

3 

6.253 

49 

62. 037 

95 

113.038 

3 

4.  644 

49 

57.  O-’C 

95 

106. 364 

4 

7.781 

50 

63.  167 

96 

114. 130 

4 

5.990 

50 

58. 164 

96 

107.425 

5 

9.237 

51 

64. 295 

97 

115.223 

5 

7.289 

51 

55  248 

97 

108.486 

6 

10. 645 

52 

65. 422 

98 

116.315 

6 

8.558 

52 

60.  331 

98 

109. 547 

7 

12.017 

53 

66.  548 

99 

117.407 

7 

9.803 

53 

61.414 

99 

110.607 

8 

13.361 

54 

67.673 

100 

118.498 

8 

11.030 

54 

62.496 

100 

111.666 

9 

14. 684 

55 

68. 796 

101 

119. 58S 

9 

12. 242 

55 

63.577 

101 

112. 726 

10 

15.987 

56 

69.918 

10 -j 

120. 678 

10 

13. 442 

56 

64. 658 

102 

113. 785 

11 

17. 275 

57 

71.  039 

3s 

121. 768 

11 

14. 631 

57 

65.  737 

103 

ii4. 844 

12 

18.  549 

58 

72. 160 

104 

122. 858 

12 

15. 812 

58 

66.816 

104 

115.903 

13 

19.812 

59 

73. 279 

105 

123. 946 

13 

16.985 

59 

67.S94 

105 

116. 961 

14 

21.064 

60 

74. 397 

106 

125.  035 

14 

18. 151 

60 

68.972 

106 

118.019 

15 

22. 307 

61 

75.514 

107 

126.  123 

15 

19.311 

61 

70.  049 

107 

119.077 

16 

23.  542 

62 

76.  630 

108 

127.211 

16 

20. 465 

62 

71.  125 

108 

120. 135 

17 

24. 769 

63 

77.  745 

109 

128. 298 

17 

21.614 

63 

72.  201 

109 

121. 192 

18 

25. 989 

64 

78. 859 

110 

129. 385 

18 

22. 759 

64 

73. 276 

110 

122. 249 

19 

27. 204 

65 

79.973 

111 

130.471 

19 

23.900 

65 

74. 350 

111 

123. 306 

20 

28.412 

66 

81.  085 

112 

131.558 

20 

25.037 

66 

75. 424 

112 

124. 363 

21 

29. 615 

67 

82. 197 

113 

132. 643 

21 

26. 171 

67 

76.498 

113 

125.419 

22 

30.813 

68 

83. 308 

114 

133.728 

22 

27.301 

68 

77.571 

114 

126.475 

23 

32.  007 

69 

84. 418 

315 

134. 813 

23 

28. 429 

69 

78.  643 

115 

127.531 

24 

33. 196 

70 

85.527 

116 

135.898 

24 

29.553 

70 

79.714 

116 

128.586 

25 

34. 382 

71 

86. 635 

117 

136. 982 

25 

30.675 

71 

80.  786 

117 

129. 642 

26 

35.563 

72 

87.  743 

118 

138. 066 

26 

31.  795 

72 

81.  856 

118 

130.697 

27 

36.  741 

73 

88.850 

119 

139. 149 

27 

32.912 

73 

82.  927 

119 

131. 751 

28 

37.916 

74 

89.956 

120 

140. 232 

28 

34. 026 

74 

83.  996 

120 

132.806 

29 

39.  087 

75 

91.061 

121 

141.315 

29 

35.  139 

75 

85.  066 

121 

133, 860 

30 

40. 256 

76 

92. 166 

122 

142. 397 

30 

36.  250 

76 

86.  134 

122 

134. 914 

31 

41.422 

77 

93.270 

123 

143. 480 

31 

37.359 

77 

87.  203 

123 

135. 968 

32 

42.  585 

78 

94. 373 

124 

144.561 

32 

38. 466 

78 

88.271 

124 

137,022 

33 

43.  745 

79 

95.476 

125 

145. 643 

33 

39.572 

79 

89.  338 

125 

138. C76 

34 

44.  903 

80 

96,578 

126 

146. 724 

34 

40.676 

80 

90.405 

126 

139.  129 

35 

46. 059 

81 

97.679 

127 

147. 805 

35 

41.  778 

81 

91.472 

127 

140,  182 

36 

47. 212 

82 

98.780 

128 

148. 885 

36 

42.879 

82 

92.538 

128 

141.  235 

37 

48.  363 

83 

99.  880 

129 

149. 965 

37 

43.978 

S3 

93.  604 

129 

142. 288 

38 

49.513 

84 

100.  980 

130 

151.045 

38 

45.076 

84 

94.  669 

130 

143. 340 

39 

50.  66C 

102. 079 

131 

152. 125 

39 

46. 173 

85 

95.  734 

131 

144. 392 

40 

5^C05 

o... 

103.  177 

132 

153.  204 

40 

47.268 

86 

96.  799 

132 

145. 444 

/,  •« 

5:...  946 

87 

104.273 

133 

154, 282 

41 

48. 363 

87 

97.863 

133 

146. 496 

42 

54.  090 

88 

lv5.  372 

134 

155.361 

42 

49.456 

88 

98.  927 

134 

147, 547 

/;3 

55.  230 

89 

106.  469 

135 

156.440 

43 

50.548 

89 

99.  990 

135 

148.599 

*+4 

56.  ''■6S 

90 

107.  565 

136 

157.517 

44 

51.639 

90 

101.  053 

136 

149. 650 

45 

91 

ICS.  660 

137 

158. 595 

45 

52.  729 

91 

102.  116 

137 

150.701 

46 

58.  640 

92 

109.  756 

138 

159.673 

46 

53.818 

92 

103.  179 

138 

151.752 
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Table  21  :  Chi-square  distribution  (Continued...) 


Left 

tail 

area  =  0, 9 

Left 

tall 

area  =  0. 8 

DF  CHI 

DP 

CHI  DF  CHI 

DF  CHI 

DF 

CHI  DF  CHI 

139  160.750 


140  161.827 

141  162.903 

142  163.979 

143  165.055 

144  166.131 

145  167.207 

146  168.282 

147  169.357 

148  170.432 

149  171.507 

150  172.581 

151  173.655 

152  174.729 

153  175.802 

154  176.875 

155  177.949 

156  179.021 

157  180.094 

158  181.167 

159  182.238 

160  183.310 

161  184.382 

162  185.453 

163  186.524 

164  187.596 

165  188.666 

166  189.737 

167  190.807 

168  191.878 

169  192.947 

170  194.017 

171  195.086 

172  196.156 

173  197.225 

174  198.294 

175  199.362 

176  200.431 

177  201.499 

178  202.567 

179  203.636 

180  204.703 

181  205.771 

182  206.838 


183  207.906 

184  208.972 

185  210.939 

186  211.  106 

187  212.172 

188  213.238 

189  214.304 

190  215.370 

191  216.436 

192  217.502 

193  218.567 

194  219.632 

195  220.698 

196  221.763 

197  222.827 

198  223.892 

199  224.956 

200  226.021 

201  227.085 

202  228.149 

203  229.213 

204  230.276 

205  231.339 

206  232.403 

207  233.466 

208  234.529 

209  235.591 

210  236.654 

211  237.717 

212  238.779 

213  23S.842 

214  240.903 

215  241.966 

216  243.027 

217  244.089 

218  245.151 

219  246.212 

220  247.273 

221  243.334 

222  249.396 

223  250.456 

224  251.517 

225  252.577 

226  253.637 


227  254.698 

228  255.758 

229  256.818 

230  257  878 

231  258.938 

232  259.997 

233  261.057 

234  262. 11? 

235  263.176 

236  264.235 

237  265.293 

238  266.353 

239  267.411 

240  268.470 

241  269.529 

242  270.587 

243  271.645 

244  272.703 

245  273.761 

246  274.819 

247  275.877 

248  276.935 

249  277.992 

250  279.050 

251  280. 107 

252  281.164 

253  282.221 

254  283.278 

255  284.335 

256  285.392 

257  286.448 

258  287.505 

259  288.562 

260  289.618 

261  290.674 

262  291.730 

263  292.786 

264  293.842 

265  294.898 

266  295.953 

267  297.010 

268  298.064 

269  299.  120 

270  300. 175 


139  152.803 

140  153.853 

141  154.904 

142  155.954 

143  157.004 

144  158.054 

145  159.103 

146  160.153 

147  161.202 

148  162.251 

149  163.300 

150  164.349 

151  165.397 

152  166.446 

153  167.494 
254  168.542 

155  169.590 

156  170.638 

157  171.686 

158  172.733 

159  173.781 

160  174.828 

161  175.875 

162  176.922 

163  177.969 

164  179.015 

165  180.062 

166  181. 108 

167  182.154 

168  183.200 

169  184.246 

170  185.292 

171  186.338 

172  187.383 

173  188.429 

174  189.474 

175  190.519 

176  191.564 

177  192.609 

178  193.654 

179  194.698 

180  195.743 

181  196.787 

182  197.832 


183  198.876 

184  199.920 

185  200.964 

186  202.007 

187  203.051 

188  204.095 

189  205. 138 
19C  206. 181 

191  207.225 

192  208.268 

193  209.311 

194  210.354 

195  211.396 

196  212.435 
19?  213.482 

198  214,524 

199  215.566 

200  216.608 

201  217.651 

202  218.693 

203  219.734 

204  220. 776 

205  221.818 

206  222.859 

207  223.901 

208  224.942 

209  225.984 

210  227.024 

211  228.066 

212  229. 107 

213  230.148 

214  231. 188 

215  232.229 

216  233.269 

217  234.310 

218  235.351 

219  236.390 

220  237.431 

221  238.472 

222  239.511 

223  240.551 

224  241.591 

225  242.631 

226  243.670 


CHI 


227  244.710 

228  245.750 

229  246.789 

230  247.328 

231  248.868 

232  249.907 

233  250.946 

234  251.985 

235  253. C24 

236  254.063 

237  255. 101 

238  256.140 

239  257. 179 

240  258.217 

241  2.59.255 

242  260.294 

243  261.332 

244  262. 371 

245  263.408 

246  264.447 

247  265.485 

248  266.522 

249  267.560 

250  268.598 

251  269.635 

252  270.673 

253  271.710 

254  272.  748 

255  273.786 

256  274.823 

257  275.860 

258  276.897 

259  277.934 

260  278.971 

261  280.008 

262  281.045 

263  282.082 

264  283. 119 

265  284. 155 

266  285. 192 

267  286.228 

268  287.264 

269  238.301 

270  289.337 
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Table  21  :  Chi-square  distribution  (Continued..  ) 


Left 

tail  area  = 

Left 

tail  area  = 

0.  8 

DF 

CHI 

DF 

CHI 

DF 

CHI 

DF 

GHI 

DF 

CHI 

DF 

CHI 

271 

301.230 

315 

347.  564 

359 

393.738 

271 

290,-  374 

315 

335.906 

359 

381.  334 

272 

302. 286 

316 

-;48.  615 

360 

394, 787 

272 

291.409 

316 

336. 940 

360 

382.366 

273 

303. 340 

317 

349.  666 

361 

395. 834 

273 

292. 445 

317 

337. 973 

361 

383. 397 

274 

304. 395 

318 

350.717 

362 

396.  882 

274 

293.481 

318 

339-007 

362 

384. 429 

275 

305. 450 

319 

351.768 

363 

397.  930 

275 

294,518 

319 

340. 039 

363 

385.459 

276 

306.505 

320 

352.819 

364 

398,978 

276 

295.553 

320 

341.073 

364 

386.491 

277 

307.559 

321 

353,869 

365 

400. 024 

277 

296.589 

321 

342. 106 

365 

387.522 

■278 

308. 614 

322 

354.919 

366 

401.071 

278 

297.625 

322 

343. 139 

366 

388.554 

279 

309. 668 

323 

355. 969 

367 

402. 119 

279 

298. 661 

323 

344. 173 

367 

389.584 

280 

310.722 

324 

357. 020 

368 

403. 167 

280 

299.697 

324 

345. 206 

368 

390.615 

281 

311.777 

325 

358.071 

369 

404, 214 

281 

300. 732 

325 

346. 239 

369 

391.  646 

282 

312.831 

326 

359. 120 

370 

405.261 

282 

301, 768 

326 

347.272 

370 

392. 678 

283 

313.  885 

327 

360. 171 

371 

406, 308 

283 

302. 803 

327 

348. 306 

371 

393. 708 

284 

314. 938 

328 

361. 222 

372 

407.356 

284 

303. 839 

.328 

349. 338 

372 

394. 740 

285 

315.992 

329 

362.271 

373 

408. 402 

285 

304. 874 

329 

350.371 

373 

395.771 

286 

317. 045 

330 

363.  322 

374 

409. 449 

286 

305.909 

330 

351.403 

374 

396.  801 

287 

318. 100 

331 

364. 371 

375 

410.496 

287 

306.944 

331 

352. 437 

375 

397. 832 

288 

319. 153 

332 

365. 421 

376 

411.543 

288 

307.980 

332 

353. 469 

376 

398. 863 

289 

320. 206 

333 

366.471 

377 

412.590 

289 

309.014 

333 

354.501 

377 

399. 894 

290 

321. 259 

334 

367.521 

37S 

413.  636 

290 

310.050 

334 

355.534 

378 

400. 924 

291 

322. 313 

335 

368.570 

379 

414. 683 

291 

311.  084 

335 

356. 567 

379 

401.  954 

292 

323. 366 

336 

369. 620 

380 

415. 729 

292 

312.119 

336 

357.599 

380 

402. 985 

293 

324. 418 

337 

370.668 

3B1 

416.776 

293 

313. 155 

337 

358. 632 

381 

404.015 

294 

325.472 

338 

371.718 

382 

417.822 

294 

314. 189 

338 

359. 664 

382 

405. 046 

295 

326. 524 

339 

3*'2.  767 

383 

418. 868 

295 

315.223 

339 

360. 697 

383 

406. 077 

296 

327.577 

340 

373.817 

384 

419.915 

296 

316.258 

340 

361.  729 

384 

407. 107 

297 

328. 630 

341 

374. 865 

385 

420. 961 

297 

317.293 

341 

362. 761 

385 

408. 138 

298 

329. 683 

342 

375.915 

386 

422. 007 

298 

318.327 

342 

363.793 

386 

409. 168 

299 

330. 735 

343 

376.964 

387 

423. 054 

299 

319,362 

343 

364. 826 

387 

410. 198 

300 

331. 788 

344 

378.012 

388 

424. 100 

300 

320.397 

344 

365. 858 

388 

411,228 

301 

332. 840 

345 

379.062 

389 

425. 146 

301 

321. 430 

345 

366.  8i,0 

389 

412. 259 

302 

333. 892 

346 

380. 110 

390 

426. 191 

302 

322. 465 

346 

367. 922 

390 

413.289 

303 

334. 945 

347 

381. 159 

391 

427. 237 

303 

323, 499 

347 

368. 954 

391 

414. 320 

304 

335.  996 

348 

382, 208 

392 

428. 284 

304 

324.533 

348 

369. 986 

392 

415.348 

305 

337. 048 

349 

383.256 

393 

429. 329 

305 

325.568 

349 

371.018 

393 

416.379 

306 

338. 100 

350 

384. 305 

394 

430.375 

306 

326. 602 

350 

372. 049 

394 

417.409 

307 

339. 152 

351 

385. 354 

395 

431. 420 

307 

327. 635 

351 

373. 082 

395 

418.439 

308 

340.  204 

352 

386. 402 

396 

432. 466 

308 

328. 670 

352 

374. 113 

396 

419.469 

309 

341.256 

353 

387.450 

397 

433.512 

309 

329. 704 

353 

375.  145 

397 

420. 499 

310 

342. 307 

354 

388.499 

398 

434. 557 

310 

330.738 

354 

376. 177 

398 

421.529 

311 

343.  358 

355 

389.547 

399 

435. 603 

311 

331.771 

355 

377.209 

399 

422.558 

312 

344. 410 

356 

390.595 

400 

436. 648 

312 

332. 805 

356 

378. 240 

400 

423. 588 

313 

345. 461 

357 

391. 642 

401 

437. 693 

313 

333. 839 

357 

379.271 

401 

424. 618 

314 

346. 512 

358 

392.690 

402 

438. 739 

314 

334. 873 

358 

380. 303 

402 

425. 648 
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